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_ CAPITOLO UNO 

Visualizzazione di una funzione 

in HI-RES con il C64 


INSERIMENTO AUTOMATICO DELLA FUNZIONE NEL LISTATO 


La funzione di cui si deve generare il grafico non può essere introdotta 
per mezzo del comando INPUT perché questo accetta solo variabili 
numeriche o di stringa e non espressioni contenenti operatori aritmetici. 
In altre parole, se tentassimo di battere 

100 INPUT "Y =" ;Y 

potremmo introdurre in Y (o in Y$) un singolo valore numerico ma mai 
una espressione contenente una variabile e segni di operazione. 

È necessario perciò usare il comando DEF FNY(X)= ..., ma sarebbe 
scomodo modificare manualmente il listato per introdurre di volta in 
volta la funzione che si vuole visualizzare sullo schermo. 

La soluzione migliore consiste in una routine che ci permetta di intro¬ 
durre automaticamente la funzione desiderata nel listato sotto la forma 
DEF FNY(X)= ... 

Si provi a lanciare la seguente routine: 

100 PRINT"(CLEAR)BATTI LA FUNZIONE 
110 INPUT"(3DOWN)Y = " ;Y$ 

120 PRINT"(CLEAR)(3DOWN)500 DEF FNY(X)="Y$ 

130 PRINT"GOTO 500(HOME)" 
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150 END 

510 PRINT'LA FUNZIONE SI TROVA ORA ALLA LINEA 500" 
Supponiamo per esempio di introdurre alla riga 110 la funzione: 
x*x-7*(x-13)+4*x-3 

sotto forma di stringa e memorizzata nella variabile Y$. 

Nelle righe 120 e 130 viene cancellato lo schermo e stampato il messag¬ 
gio seguente: 


2 RIGHE VUOTE j 

POSIZIONE-► 

CURSORE 


II2A2Y. 

500 22F HIY(X)=»*X-7*(x-I3)+4*X-3 
GOTO 500 


Fig. 1.1 


Infatti il computer stampa i messaggi contenuti in 120 e 130 (rispetti¬ 
vamente sulla quarta e quinta riga di schermo), poi il cursore viene 
riportato in HOME e viene stampato il messaggio READY di fine pro¬ 
gramma. 

Se a questo punto battete per due volte consecutive il tasto RETURN 
vengono rese operative le scritte che passano sotto il cursore. 

Quindi viene prima introdotta la riga 500 nel listato (o viene modificata 
nel caso in cui fosse già presente) e successivamente si riprende a 
scorrere dalla riga 500. 

È facile verificare, con un LIST, che effettivamente la riga 500 è stata 
introdotta nel programma in modo definitivo. 

Questo procedimento può essere migliorato forzando tre RETURN 
automatici ed eliminando così la necessità di eseguire manualmente tale 
operazione. 

Il computer ha la capacità di memorizzare l’ordine con cui vengono 
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premuti più tasti in successione. La porzione di memoria adibita a tale 
impiego si chiama buffer di tastiera e inizia alla locazione 631. 

Nella locazione 198 viene invece precisato quanti caratteri il computer 
deve prelevare dal buffer di tastiera e rendere operanti (fino ad un 
massimo di 10 caratteri). 

A questo punto possiamo sfruttare il buffer di tastiera aggiungendo al 
listato precedente la riga 

150 FOR K=0 TO 2:POKE 631 +K,13:NEXT:POKE 198,3 

Grazie al ciclo FOR...NEXT vengono introdotti tre RETURN (CHR- 
$(13)=RETURN) nel buffer di tastiera e quindi viene precisato al 
calcolatore di prelevare tre caratteri dal buffer (infatti nella locazione 
198 viene detto al calcolatore quanti caratteri deve prelevare dal buffer 
di tastiera). Abbiamo aggiunto un RETURN più del necessario per 
maggior sicurezza e per rendere operativa una eventuale altra riga. 
Provate a far girare la routine e vi accorgerete di quanto essa sia efficace 
e suscettibile di applicazioni anche per altri impieghi. 

È buona norma riportare la locazione 198 in condizioni normali dopo 
aver sfruttato queste routine, con il comando 

POKE 198,0 


Finalmente siamo in grado di realizzare una routine definitiva che ci 
permette di introdurre nel listato la funzione da presentare in grafico, di 
stabilire gli intervalli di visualizzazione e di precisare la presenza di 
eventuali asintoti verticali. 

La necessità di segnalare al computer la presenza degli asintoti nasce dal 
fatto che nelle ascisse corrispondenti l’ordinata è infinita e perciò nel 
calcolo dell’ordinata stessa il calcolatore incorrerebbe nell’errore di 
divisione per zero o nell’errore di overflow. 

Conoscendo invece tali ascisse eviteremo il calcolo dell’ordinata facen¬ 
do saltare il programma al tracciamento dell’asintoto. 

100 CLR :PRINT"C CLEARIC BLEUI":POKE 53281,15 : POKE 53280,8 
110 PRINT"CCLEARIBATTI LR FUNZIONE " 

120 INPUT "C3 DONNIV « ",V* 

138 PRINT"CCLEARIC3 DONNI500 DEF FNY(X)*"V# 

140 PRINT"GOTO 500CHOMEI" 

150 FOR K*0 TO 2 POKE 631 *K ,13 : NEXT'POKE 198,3 
160 END 

505 POKE 198,0 : PRINT"CCLERR]OK - ORA ") 

510 PRINT"SCEGLI OLI INTERVALLI DI" 

520 PRINT TRB(5)"VISUALIZZAZIONE :" 

530 INPUT "[2 DONNI* (SINISTRA)»",XI 
540 INPUT "X (DESTRA)*",X2 
550 IF X1>=X2 THEN 510 
560 INPUT "CDOWNIV (BASSA)*".;VI 
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570 INPUT "V <RLTR)*";V2 
580 IF Y1>*Y2 THEN 510 

590 PRINT"CB0WN]CI SONO ASINTOTI VERTICALI ? 

600 GOSUB 15080 PRINT GSIF Q*="N" THEN 660 
610 PRINT"[DOWN3QUANTI ? ", 

620 GOSUB 15000:PRINT Q*:Nfi*Q 
630 PRINT : FOR K*1 TO NA 

640 PRINT"ASCISSA"K"* ASINTOTO : INPUT ASCIO 
650 NEXT K 
660 END 

15000 REM - 

15010 REM ROUTINE DI ATTESA 

15020 REM - 

15030 GET Q$ IF Q$="" THEN 15030 
15040 Q=VALCQ*) RETURN 


Non abbiamo ritenuto necessario dimensionare la variabile AS(K), 
perché consideriamo assolutamente improbabile che sia necessario vi¬ 
sualizzare più di 10 asintoti verticali contemporaneamente e, come è 
noto, nei computer Commodore non è necessario dimensionare le 
variabili con indice quando il dimensionamento non supera il valore 
dieci. 

La subroutine 15000 serve per creare un ciclo che ... tenga impegnato il 
computer in attesa che venga premuto un tasto. 


ALTA RISOLUZIONE E TRACCIAMENTO 


Ogni carattere che appare sullo schermo è costituito da una matrice di 8 
x 8 punti (o pixel) opportunamente accesi o spenti (cioè con le stesse 
dimensioni del cursore lampeggiante). 

In condizioni normali lo schermo contiene mille caratteri e quindi è 
strutturato con: 

8 x 8 x 1000 = 64000 punti 


Si può accedere all’alta risoluzione, cioè alla gestione dei singoli pixel, 
per mezzo della seguente routine: 


100 POKE 53272,PEEK(53272) OR 8 
110 POKE 53265,PEEK(53265) OR 32 
120 FORK=8192TO 8192+7999 
130 POKE K,0 
140 NEXT K 

150 FOR K=1024 TO 2023 
160 POKE K,3 
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170 NEXTK 
180 RETURN 

Nella riga 100 viene forzato un 1 nel quarto bit della locazione 53272 
(registro interno del vìe che controlla l’inizio della matrice video). 

In 110 viene forzato un 1 nel sesto bit della locazione 53265 (registro 
interno del vie adibito a interruttore per attivare l’alta risoluzione). 
Nel primo ciclo FOR...NEXT viene pulito tutto l’intervallo di memoria 
riservato alla memoria di schermo in alta risoluzione. Essa inizia alla 
locazione 8192 e termina dopo 8000 byte, in quanto ogni locazione 
raccoglie i dati di otto pixel (1 pixel per ciascun bit). 

Infine con il secondo ciclo FOR...NEXT coloriamo in cyan lo sfondo 
dello schermo. Infatti in alta risoluzione l’intervallo di memoria che 
comunemente è riservato ai caratteri di schermo viene trasformato in 
memoria del colore. 

Si noti come, per questa ragione, non è possibile inserire più di un colore 
all’interno di un quadratino di 8 x 8 pixel. 

Questa routine ha lo svantaggio di essere piuttosto lenta ed è quindi 
opportuno sostituirla con una identica subroutine in linguaggio mac¬ 
china. 


20000 REM - 

20010 REM ROUTINE PER RTTIVRZ. HI-RES 

20020 REM - 

20030 F0R K> 49152 T0 49222 

20040 READ N 

20050 P0KE K,N 

20070 NEXT K 

20080 RETURN 

20090 DATA 173,24/208,9,6.141,24,208,173,17,208,9,32,141,17,208 
.169,0,133 

20100 IiATfl 251,169,32,133,252,160,0,169,0,145,251,200,192,0,208 
249,230 

20110 URTA 252,169,64,197,252,208,239,169,0,133,251,169,4,133,2 
52,160,0 

20120 URTA 169,3,145,251,200,192,0,288,249,230,252,169,8,197,25 
2,208,239.96 

Con l’attivazione dell’alta risoluzione possiamo ora operare su un reti¬ 
colato di pixel ordinati nel modo seguente. 
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Fig. 1.2 


Il byte zero è quello corrispondente alla locazione 8192. Si noti come il 
primo pixel che si trova nell’angolo superiore sinistro dello schermo 
corrisponde al bit più alto del byte zero mentre l’ultimo pixel in basso a 
destra corrisponde al bit più basso della locazione 8192 + 7999. 

Le coordinate di ciascun pixel sono fissate con riferimento all’angolo 
superiore sinistro dello schermo. Perciò il pixel di coordinate 100,50 
corrisponde al punto che si ottiene partendo dall’angolo suddetto e 
spostandoci verso destra di cento pixel e verso il basso di 50. 
Osservando la figura risulta evidente la disposizione apparentemente 
innaturale con cui procedono i pixel mano a mano che le ascisse o le 
ordinate aumentano. È necessario perciò strutturare comandi che faci¬ 
litino l’individuazione del pixel da accendere. 

Supponiamo di voler accendere un unico pixel di coordinate x,y. 

INT (Y/8) Serve per determinare a quale delle 25 righe 

appartiene il pixel. 

INT (X/8) Serve per determinare a quale delle 40 colonne 

appartiene il pixel. 

Y AND 7 Serve per determinare a quale degli otto byte 

(della matrice quadrata di 8 x 8 bit) appartiene il 
pixel. 
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7—(X AND 7) Serve per determinare a quale bit del byte appar¬ 

tiene il pixel. 

Per l’accensione del pixel desiderato possiamo quindi usare la subrou¬ 
tine seguente: 

25030 IF X<0 OR X>319 OR Y<0 OR Y>199 THEN 25060 
25040 BY=8192*320*INT(Y/8)+8*!NT(X/8)+(Y AND 7) 

25050 POKE BY,PEEK(BY) OR 2*(7-(X AND 7)) 

25060 RETURN 

In 25030 viene effettuato un controllo per evitare l’uso di coordinate 
illegali. 

Nella riga 25040 si determina in quale byte è contenuto il bit corrispon¬ 
dente al pixel da accendere. 

Infine nella riga 25050 viene forzato un 1 in corrispondenza del bit 
desiderato. 

Per disattivare l’alta risoluzione e tornare in bassa risoluzione, usate i 
due comandi seguenti: 

30030 POKE 53272,PEEK(53272) AND 247 
30040 POKE 53265,PEEK(53265) AND 223 

dove 247 e 223 sono rispettivamente i complementi di 8 e 32 che con 
l’operazione logica AND forzano uno zero nei due registri già visti nella 
fase di attivazione dell’alta risoluzione e che erano stati portati a 1. 
Abbiamo così ristabilito le condizioni iniziali e siamo tornati in bassa 
risoluzione. 

Il listato precedente può quindi essere integrato nel modo seguente 


100 G0SUB 20000 : X*160 : V«100 : SVS49152 G0SUB 25000 
110 GET Q*IF G*="" THEN 110 
120 G0SUB 30000 END 

20000 REM - 

20010 REM ROUTINE PER RTTIVHZ. HI-RES 

20020 REM - 

20030 F0R K= 49152 T0 49222 

20040 RERD N 

20050 POKE K,N 

20070 NEXT K 

20080 RETURN 

20090 DHTfl 173,24,208,9,8,141,24,208,173,17,208,9,32,141,17,208 
■ 169,0•133 

20100 DhTFi 251,169,32,133,252,160,0,169,0,145,251,200,192,0,208 
■249,230 

20110 DHTFI 252,169,64,197,252,208,239,169,0,133,251,169,4,133,2 
52,160,0 

20120 DATA 169,3,145,251,200,192,0,208,249,230,252,169,8,197,25 
2,208,239,96 
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25000 REM - 

25010 REM ROUTINE PER COMANDO PLOT 

25020 REM -^- 

25030 IF X<0 OR >0319 OR V<0 OR V>199 THEN 25060 
25040 BV='8192+320#INT<V/8>+8i(fINT<XZ8) + <V AND 7) 

25050 POKE BV,PEEK<BV) OR 2t<7-(X AND 7)) 

25060 RETURN 

30000 REM - 

30010 REM ROUTINE PER DISATTIV. HI-RES 

30020 REM - 

30030 POKE 53272,PEEK<53272) AND 247 
30040 POKE 53265,PEEK<53265) AND 223 
30050 PRINT"CCLEARl" RETURN 


Le righe 100, 110 e 120 servono solo per controllare il buon funziona¬ 
mento delle subroutine. Se tutto è stato digitato correttamente, si dovrà 
accendere un unico pixel al centro dello schermo (coordinate 160,100). 


ELABORAZIONE E TRACCIAMENTO PUNTI 


Con l’alta risoluzione lo schermo viene trasformato in una griglia di 
64000 punti disposti su 200 righe e 320 colonne. Ciascuno di questi punti 
viene chiamato pixel e il pixel di coordinate 0,0 corrisponde all’angolo 
superiore sinistro dello schermo. L’ultimo pixel è invece quello con 
coordinate 319, 199 e corrisponde all’angolo inferiore destro. 

Ci occuperemo di funzioni (algebriche o trascendenti) del tipo: 

y = /(*) 

cioè funzioni univoche con al massimo un solo valore per ogni ascissa. 
Il grafico della funzione si otterrà con un ciclo FOR...NEXT per le 
ascisse che vanno da zero a 319. Per ciascuna di queste ascisse sarà 
calcolata la corrispondente ordinata e poi questa verrà visualizzata con 
l’accensione del relativo pixel. 

Un primo problema importante da chiarire è la differenza tra coordina¬ 
te assolute e relative. Infatti noi scegliamo un opportuno intervallo di 
visualizzazione sia per le ascisse che per le ordinate, aprendo così 
un’ipotetica “finestra” entro la quale avremo la visualizzazione della 
funzione. Le coordinate relative a tale finestra saranno appunto le 
coordinate relative e generalmente queste saranno differenti dalle coor¬ 
dinate effettive o assolute della funzione. 

Indicheremo con x u x 2 , y t , y 2 gli estremi della finestra, come indicato in 
figura. 
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Per esempio il punto A della parabola ha coordinate assolute 1,2 e 
coordinate relative 160,100. 

Il passaggio da coordinate assolute (x 0 , y 0 ) a coordinate relative (x, y) si 
ottiene con le formule seguenti: 


x = (x,) - X\)IDx 

y = Cvo - y\ )/Dy 

dove x e y rappresentano rispettivamente il numero di unità assolute che 
intercorrono fra un pixel e il successivo 

Dx = (x 2 — Ar!)/320 
Dy = (y 2 - yi)/200 

Le formule inverse per passare da coordinate relative a coordinate 
assolute sono invece: 

jc 0 = + xDx 

)>0 = y\ + yDy 

Una volta fissati gli intervalli di visualizzazione, potremo quindi calco¬ 
lare le coordinate relative da riportare in grafico con la routine se¬ 
guente. 
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1030 SYS 49152 : REM ATTIVAZIONE HI-RES 
1040 DX = (X2-X1)/320 : DY = (Y1-Y2)/200 
1060 FOR X=0 TO 319 
1070 X0—XI +X*DX 
1080 Y0=FNY(X0) 

1090 IF Y0<=Y1 THEN 1130 
1100 IFY0>=Y2 THEN 1130 
1110 Y=200-(Y0-Y1)/DY 

1120 GOSUB 25000: REM SUBROUTINE DI PLOTTAMENTO 
1130 NEXT X 

1140 GOSUB 15000: REM SUBROUTINE DI ATTESA 
1150 GOSUB 30000: REM SUBROUTINE DI RITORNO IN BASSA 
RISOLUZIONE 

Nella riga 1110 viene prima calcolata Tordinata relativa ((Y0-Y1 )/DY) e 
quindi viene determinato il complemento a 200 perché nelle coordinate 
assolute la Y aumenta andando verso l’alto, mentre in quelle relative la 
Y aumenta andando verso il basso. 

Le righe 1090 e 1100 servono a saltare la routine di tracciamento nei 
casi in cui l’ordinata non risulti dentro la finestra di visualizzazione. 


LISTATO DEFINITIVO 


Possiamo ora procedere alla elaborazione del listato finale. 

Nella parte introduttiva (fino alla riga 650) viene introdotta la funzione, 
vengono scelti gli intervalli di visualizzazione e viene notificata al com¬ 
puter l’ascissa relativa agli asintoti verticali. 

Si passa alla routine di elaborazione con una SYS di attivazione dell’alta 
risoluzione, e con la determinazione dei parametri X e Y. 

Nella riga 1050 si ha l’invio al tracciamento degli assi cartesiani. 
Nelle righe comprese fra la 1060 e la 1130, con un ciclo FOR...NEXT si 
prendono in considerazione tutte le ascisse relative comprese tra 0 e 
319; per ciascuna viene determinato il corrispondente valore assoluto 
X0 (riga 1070) e si ha il passaggio per la subroutine 13000 in cui si ha un 
test per controllare se tale ascissa assoluta coincide con uno degli 
asintoti. In caso affermativo si passa alla riga 13060 per il tracciamento 
dell’asintoto stesso, l’incremento dell’ascissa relativa, il calcolo della 
nuova ascissa assoluta e il ritorno dalla subroutine (riga 13100). 

In caso negativo si ha il ritorno immediato nella linea 13050. 
Ritornando nel programma principale in 1080, si ha il calcolo dell’ordi¬ 
nata assoluta mentre, nelle due linee successive viene controllato se tale 
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ordinata cade nella finestra della visualizzazione. In caso negativo si ha 
il salto alla riga 1130. In caso affermativo nella riga 1110 viene calcolata 
l’ordinata relativa, più precisamente il suo complemento a 200. Si passa 
quindi alla subroutine 25000 di tracciamento. 

Infine in 15000 si ha la routine di attesa e in 30000 la disattivazione 
dell’alta risoluzione. 


100 CLR : PPINT"[CLEAR3 CBLEU3" : P0KE 53281,15 :P0KE 53280,8 G0S 
IIB 20000 

110 PRINT"[CLEAR3BATTI LA FUNZIONE " 

120 INPUT "C3 DOWN3V = ";V* 

130 PRINT"CCLEARK3 DOWN3500 DEF FNV(X)*"V* 

140 PRINT"GOTO 500CHOME3” 

150 F0R K=0 T0 2 : P0KE 631+K,13 :NEXT:P0KE 198,3 
160 END 

500 DEF FNV(X > = < X-2 V < X+2) 

505 POKE 198,0 PRINT"COLEAR3OK - ORA "; 

510 PRINT"SCEGLI GLI INTERVALLI DI" 

520 PRINT TAB<5)"VISUALIZZAZIONE " 

530 INPUT "[2 DOWN3X (SINISTRA)*";XI 

540 INPUT "X (DESTRA)*";X2 

550 IF X1>*X2 THEH 510 

560 INPUT "CD0WN3V (BASSA)*";VI 

570 INPUT "V (ALTA)*";V2 

580 IF V1>*V2 THEN 510 

590 PRINT"CD0WN3CI SONO ASINTOTI VERTICALI > ", 

6Ù0 00SUB 15000 PRINT 0*1F G*="N" THEN 1000 
610 PRINT"[DOWN3GUANTI '' ": 

620 GOSUB 15000 PRINT G*:NA=G 
630 PRINT : FOR K*1 TO Nfl 

640 PRINT"ASCISSA"K"* ASINTOTO =", INPUT ASCIO 
650 NEXT K 

1000 REM - 

1010 REM ELABORAZIONE 

1020 REM - 

1030 SVS 49152 

1040 DX=CX2-X1)X320 DV*CV2-V1>2200 
1050 GOSUB 10000 GOSUB 11000 
1060 FOR X*0- TO 319 
1870 X0“X1+X*DX : GOSUB 13080 
1880 V0*FNVCX0) 

1090 IF V0<=V1 THEN 1130 

1180 IF V0>*V2 THEN 1130 

1110 V=200-(V0-V1VDV 

1120 GOSUB 25000 

1130 NEXT X 

1140 GOSUB 15000 

1150 GOSUB 38000 

2808 REM - 

2010 REM MESSAGGI FINALI 

2820 REM - 

2030 PRINT"VUOI CAMBIARE GLI INTERVALLI DI" 

2040 PRINT"VISUALIZZAZIONE ?" 

2050 GOSUB 15000 

2060 IF 0$*"S" THEN PRINT"CCLEAR3 ";:G0T0 510 

2070 PRINT"[DOWN3VUOI CAMBIARE LA FUNZIONE ?" 

2080 GOSUB 15000 

2090 IF G**"S" THEN RUN 

2100 PRINT"CD0WN30K - FINE LAVORO." 

2110 END 

10000 REM - 
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10010 REM ASSE X 

10020 REM - 

10030 IF SGN<V1)■ S0N<V2) THEN RETURN 
10040 V« I NT < V2XDV+. 5 ) 

10050 FOR X=0 TO 319 STEP 6 
10060 GOSUB 25000 
10070 NEXT X RETURN 

11000 REM - 

11010 REM ASSE V 

11020 REM -*- 

11030 IF SGN<X1>* SGNCX2) THEN RETURN 
11040 X=-1NT < X1/DX+.5 > 

11050 FOR V*0 TO 199 STEP 6 
11060 GOSUB 25000 
11070 NEXT V RETURN 

13000 REM - 

13010 REM ASINTOTI 

13020 REM - 

13030 FOR J=1 TO Nfi 

13040 IF xe-ftS<J> THEN 13060 

13050 NEXT J RETURN 

13060 FOR V=0 TO 199 STEP 3 

13070 GOSUB 25000 

13080 NEXT V 

13090 X=X+1 X0=X1+X*DX 

13100 RETURN 

15000 REM - 

15010 REM ROUTINE DI ATTESA 

15020 REM - 

15030 GET 0* IF Q$="" THEN 15030 
15040 Q=VAL (OS) RETURN 

20000 REM - 

20010 REM ROUTINE PER ATTIVAZ. HI-RES 

20020 REM - 

20030 FOR K> 49152 TO 49222 

20040 READ N 

20050 POKE K,N 

20070 NEXT K 

20080 RETURN 

28090 URTA 173,24,208,9,8,141,24,208,173,17,208,9,32,141,17,208 
,169,8,133 

20100 DATA 251,169,32,133,252,160,0,169,0,145,251,200,192,0,208 
,249,230 

20110 DATA 252,169,64,197,252,208,239,169,0,133,251,169,4,133,2 
52,160,0 

20120 DATA 169,3,145,251,200,192,0,208,249,230,252,169,8,197,25 
2,208,239,96 

25000 REM - 

25010 REM ROUTINE PER COMANDO PLOT 

25020 REM - 

25030 IF X<0 OR X>319 OR V<0 OR V>199 THEN 25060 
25048 BV=8192+328*1NT<V/8)+8*1NT<X/8)+(V AND 7) 

25050 POKE BV,PEEK<BV) OR 2T(7-<X AND 7)) 

25060 RETURN 

38080 REM - 

38010 REM ROUTINE PER DISATTIV. HI-RES 

30020 REM - 

30030 POKE 53272,PEEK<53272) AND 247 
30040 POKE 53265,PEEK(53265) AND 223 
30050 PRINT"[OLEARI" RETURN 


Per il tracciamento degli assi coordinati nelle apposite subroutine 10000 
e 11000 si ha prima il controllo che i segni dei due estremi di visualizza- 
















VISUALIZZAZIONE DI UNA FUNZIONE IN HI-RES CON IL C64 


zione sono differenti: infatti se essi sono concordi l’asse opposto cade 
sicuramente fuori della finestra e non deve essere visualiz z ato. 

In tal caso si ha il RETURN immediato mentre, in caso contrario, viene 
calcolata la seconda coordinata relativa e con un ciclo FOR...NEXT e 
l’uso della subroutine 25000 viene realizzato il tracciamento dell’asse in 
questione. 

Il programma principale termina con le righe da 2000 a 2110 contenenti 
le richieste finali. 

In caso si vogliano cambiare gli intervalli di visualizzazione viene evitato 
il passaggio di introduzione della funzione e di lettura dei DATA occor¬ 
renti per la routine dell’alta risoluzione. 



CAPITOLO DUE 

Stampa di ima funzione 
con il C64 


ROUTINE DI HARD COPY 


Le funzioni visualizzate sullo schermo possono essere trasferite su carta 
con una opportuna routine in linguaggio macchina, sistemata in memo¬ 
ria a partire dalla locazione 52224 fino alla locazione 52667, e che 
funziona con la stampante MPS-801 della Commodore. 

La stampa viene eseguita trasversalmente nello stesso verso di scorri¬ 
mento del foglio, ed è possibile ottenere una dilatazione o meno lungo 
l’asse x con i comandi: 


POKE 2,2 stampa amplificata 
POKE 2,1 stampa ristretta 

100 CLP ■ F'RINT " [CLEAR ] CBLEU3 " '■ POKE 53281,15 : POKE 53230,8 

105 PRINTTAB(Z)"C2 DOWN3ATTENDERE PER FAVORE.":GOSUB 200 

00GOSUB 32100 

110 PRINT"[CLEAR3BATTI LA FUNZIONE : " 

120 INPUT "[3 DOWNIV = ",V$ 

130 PRINT"CCLEARI[3 DOWN3500 DEF FNV(X)="V$ 

140 PRINT"00T0 500[HOME3" 

150 FOR K=0 TO 2 POKE 631+K,13 : NEXT : POKE 138,3 
160 END 

500 DEF FNV<X>*2#SIN<X)+< COS <X)12) 

505 POKE 198,0 PRINT"tCLEAR3OK - ORA ", 

510 PRINT-SCEOLI GLI INTERVALLI DI" 

520 PRINT TAB(5)"VISUALIZZAZIONE 
530 INPUT "12 DOWN3X (SINISTRA XI 
540 INPUT "X (DESTRA>=",X2 
550 IF XI >=X2 THEH 510 
560 INPUT "CDOWN3V (BASSA)*";VI 
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570 INPUT "V ■ALTA>=",V2 
580 IF VI>=V2 THEH 510 

580 PR1NT"CD0NN1CI SONO ASINTOTI VERTICALI " " • 

6.00 GOSUB 15000 F'RINTQJ • IF 0*="N" THEH NA=-1 GOTO 1000 
610 PRINT"[DONNIQUANTI ? "j 
620 GOSUB 15000 PPINTUÌ NA=Q 
630 PPINT : FOP K=1 TO HA 

640 F'PINT"ASC18SA"K"* ASINTOTO INPUT ASCIO 

650 ne: :t k 

1000 REM - 

1010 REM ELABORAZIONE 

1020 REM - 

1030 SYS 4SI52 

1040 DX= «:X2-X1 >/320 DV= < V2-V1 ) /200 
1050 GOSUB 10000 GOSUB 11000 
1060 FOR X=0 TO 319 
1070 X8=X1+X*DX:GOSUB 13000 
1080 V0=FNV C X0 > 

1090 IF V0OV1 THEH 1130 

1100 IF V0>=V2 THEN 1130 

1110 V=200-<V0-V1>/DV 
1120 GOSUB 25000 
1130 NEXT X 
1140 GOSUB 15000 

1145 IF THEH GOSUB 32800 

1150 GOSUB 30000 

2000 REM - 

2018 REM MESSAGGI FINALI 

2020 REM - 

20SA PRINT"VUOI CAMBIARE GLI INTERVALLI DI" 

2040 PRINT"VISUALIZZAZIONE ?“ 

2050 GOSl IB 15000 

2060 IF QI="S" THEN PRINT"CCLEAR] "OG0T0 510 

2070 PRINT"[DONNIVUOI CAMBIARE LA FUNZIONE 9 " 

2080 GOSUB 15000 

2090 ÌF QI*"S" THEN RUN 

2100 PRINT"[DONNIOK - FINE LAVORO." 

2110 END 

10000 REM - 

18010 REM ASSE X 

10020 REM - 

10030 IF SGN<V1)= SGNCV2) THEN RETURN 

10040 V=INT CV2/DV+.5 > 

10050 FOR X=0 TO 319 STEP 6 
10060 GOSUB 25000 
10070 NEXT X : RETURN 

11000 REM - 

11010 REM ASSE V 

11020 REM - 

11030 IF SGN<XI >= SGNCX2 ') THEN RETURN 
11040 X=-1NT < X1/DX+.5 > 

11050 FOR V=0 TO 199 STEP 6 
1106fi GOSUB 25000 
1107O NEXT V RETURN 

13800 REM - 

13010 REM ASINTOTI 

13020 REM - 

13025 IF NA=-1 THEN 13100 

13030 FOR J=1 TO NA 

13040 IF X0*AS<J> THEN 13068 

13050 NEXT J RETURN 

13060 FOR V=0 TO 199 STEP 3 

13*78 GOSUB 25080 

13*80 NEXT V 

13090 X=X+1 : XO=X1+X4DX 
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13100 RETIJRH 

15000 REM - 

15010 REM ROUTINE DI ATTESA 

15020 REM - 

15030 GET Ùf IF Df="" THEN 15030 
15040 G=VAL (.0$ > : RETURN 

20000 REM - 

20010 REM ROUTINE PER ATTIVAZ. HI-RES 

20020 REM - 

20030 POR ¥- 49152 TO 49222 

20040 READ N 

20050 POKE K , N : ND=ND+N 

20070 NEXT K 

20080 RETURN 

20090 DATA 173,24,208,9,8,141,24,208,173,17,208,9,32,141,17,208 
,169,0,133 

20100 DATA 251,169,32,133,252,160,0,169,0,145.251,200,192,0.208 
,249,230 

20110 DATA 252,169,64,197,252,208,239,169,0,133,251,169,4,133,2 
52,160,0 

20120 DATA 169,3,145,251,200,192,0,208,249,230,252,169,8,197,25 
2,208,239,96 

25000 REM - 

25010 REM ROUTINE PER COMANDO PLOT 

25020 REM - 

25030 IF XC0 OR X>319 OR Y<0 OR V>199 THEN 25060 
25040 BV=8192+320#INTC'tV8)+8*INT<X/8) + ( V AND 7) 

25050 POKE BV,PEEKCBV) OR 2TC7-CX AND 7)> 

25060 RETURN 

30000 REM - 

30010 REM ROUTINE PER DISATTIV. HI-RES 

30020 REM - 

30030 POKE 53272,PEEKC53272) AND 247 
30040 POKE 53265,PEEKC53265) AND 223 
30050 PRINT”CCLERR]" RETURN 
32000 REM - 


32010 REM STAMPA 

32020 REM - 

32030 FOR X=A TO 319 STEP 2 : Y=0■GOSUB 25000 V»199:OOSUB 

5000 NEXT X 

32040 FOR^ V=0^ TO 199 STEP 2X=0GOSUB 25000 X*319 OOSUB 

32050 SVS 52224 
32060 RETURN 

32070 REM - 


32080 REM ROUTINE PER STAMPA 
32090 REM - 

32100 POKE 2,2 REM PER STAMPA RISTRETTA SOSTITUIRE CON POKE 

32110 FOR K=52224 TO 52677 

32120 READ N:POKE K,N ND=ND+N NEXT K 

32130 IF NDO70923 THEN PRINT"tDOWN]ERRORE NEI DATA" END 
32140 RETURN 

32270 DATA 173,0,221,41,3,73 
32280 DATA 3,160,6,10,136,208 
32290 DATA 252,141,132,3,173,24 
32300 DATA 208,41,8,240,9.24 


32310 DATA 
32320 DATA 
32330 DATA 
32340 DATA 
32350 DATA 
32360 DATA 
32370 DATA 
32380 DATA 


169,32,109,132,3, 141 
132,3,169,0,32,189 

255.169.4.170.160.255 

32.186.255.32.192.255 
162,4.32,201,255,176 
3,76,61,204,76,32 

205.169.8.32.210.255 
169,13,32,210,255,162 


2 

2 
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32398 DATA 
32488 DATA 
32418 DATA 
32420 DATA 
32430 DATA 
32440 DATA 
32450 DATA 
32468 DATA 
32470 URTA 
32488 DATA 
32490 DATA 
32500 DATA 
32510 DATA 
32520 DATA 
32530 DATA 
32540 DATA 
32550 DATA 
32560 DATA 
32578 DATA 
32580 DATA 
32590 DATA 
32680 DATA 
32610 DATA 
32620 DATA 
32630 DATA 
32640 DATA 
32650 DATA 
32660 DATA 
32670 DATA 
32680 DATA 
32690 DATA 
32700 DATA 
32710 DATA 
32720 DATA 
32730 DATA 
32740 DATA 
32758 DATA 
32760 DATA 
32778 DATA 
32780 DATA 
32790 DATA 
32800 DATA 
32810 DATA 
32820 DATA 
32830 DATA 
32840 DATA 
32850 DATA 
32860 DATA 
32870 DATA 
32880 DATA 
32890 DATA 
32900 DATA 
32910 DATA 
32920 DATA 
32930 DATA 
32940 DATA 
32950 DATA 
32960 DATA 
32970 DATA 
32980 DATA 
32990 DATA 
33000 DATA 
33010 DATA 
33020 DATA 


0,169,1,141,198,205 
169,0,141,199,205,169 
R,141,200,205,169,199 

141.201.205.32.225.255 
208,3,76,32,205,138 

72.152.72.32.41.205 

104.168.104.170.173.205 
205, 45,203,205,240,12 

173.202.205.13.198.205 

141.282.205.76.143.204 

173.198.205.73.255.45 
202,205,141,202,285,14 
198,205,173,198,205,201 

128.240.20.24.173.199 

205.105.1.141.199.205 
173,200,205,105,0,141 

200.205.76.93.204.173 

202.205.9.128.224.45 

144.10.173.202.205.41 

31.9.128.141.202.205 
142,207,205,166,2,142 

206.205.168.32.210.255 

152,286,206,205.208,246 

174.207.205.169.1.141 
198,205,169,0,141,202 
205,56,173,199,205,233 

6.141.199.205.173.200 
205,233,0,141,200,205 

206.281.205.173.201.205 
201,255,240,3,76,93 
204,224,45,176,31,24 

173.199.205.105.7.141 
199,205,173,200,205,105 
0,141,200,205,232.169 
199-141,201,205,169,13 

32.210.255.76.93.204 
169,13,32,210,255,32 

231.255.96.173.201.205 

41.7.141.204.205.173 
201,205,74,74,74,168 
185,146,205,133,251,185 
172.205,133,252,24,165 
251,109,204,205,133,251 
165,252,105.0,133.252 
24,173,132,3,101.252 

133.252.173.199.205.41 
7 . 73 , 7 , 168 ,200,169 
0,56,42,136,208,252 

141.203.205.24.173.200 
205.101,252, 1 33,252, 173 
199,205,41,248,168, 138 

72.120.162.52.134.1 

177.251.162.55.134.1 

88.168.104.170.152.45 
203,205,141,205,205,96 
0.64.128.192,0,64 
128,192,0,64,128,192 
0,64,128,192,0,64 
128,192,0,64.128, 192 
0,64,0,1.2,3 
5,6,7,8,10,11 

12,13,15,16,17,18 
20 ,21,22', 23,25,26 
27,28,30,31 
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Dopo il tracciamento della funzione sullo schermo (riga 1140) si ha 
l’invio nella routine d’attesa e nel caso che venga premuto il tasto <— 
viene eseguito un salto alla subroutine di riga 32000 per l’esecuzione 
della stampa. 

Nelle righe 32030 - 32040 viene tracciata una cornice, il cui scopo è 
delimitare sulla carta la finestra di visualizzazione della funzione. 

In 32050 viene attivata la routine di hard copy, e infine si ha il RETURN 
per i messaggi finali. 

Si noti che talvolta il programma sembra fermarsi inspiegabilmente per 
alcuni secondi quando l’ascissa assume il valore zero; ciò è dovuto a una 
confusione di variabili che si verifica quando non ci sono asintoti verti¬ 
cali. 

Questo inconveniente può essere eliminato ponendo NA = — 1 (riga 
600) nel caso in cui non vi siano asintoti e aggiungendo la riga 13025 che 
realizza un salto impedendo il tracciamento dell’asintoto in corrispon¬ 
denza dell’asse y. 

Talvolta gli asintoti verticali possono non essere tracciati, e questo 
avviene quando a partire dall’ascissa minima jc t con incrementi succes¬ 
sivi àjc l’asintoto stesso viene a trovarsi in una posizione intermedia. 
Un aspetto importante da sottolineare è l’impossibilità (per esempio y 
= tg (x); Y = sec (X); Y = X/(X f 2 - 3)...) di segnalare per mezzo 
dell’INPUT asintoti verticali aventi per ascissa numeri irrazionali (X = 3t 
/ 2; X = SQR(3)...). 

Per fortuna non è un problema: partendo dall’ascissa forzatamente 
razionale X\ e con incrementi anche essi razionali, non si potrà mai 
pervenire all’elaborazione di una ascissa irrazionale. Quindi in ogni 
caso tali asintoti verticali verranno ignorati senza che sia necessario 
notificarli al calcolatore prima del tracciamento. 


AVVERTENZE PER LA STAMPA 


Il programma traccia la funzione con la massima precisione possibile 
senza tralasciare alcun punto e quindi con una velocità che talvolta può 
essere desiderabile aumentare. 

La cornice eseguita all’inizio della fase di stampa ha lo scopo di delimi¬ 
tare il campo di tracciamento su carta con un effetto che riteniamo 
esteticamente valido. 

È però possibile velocizzare il tracciamento del grafico eliminando 
completamente tale cornice, sopprimendo le righe 32030 - 32040, 
oppure sostituendole con: 
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32030 FOR X=0 TO 319 STEP 6:Y=0:GOSUB 
25000:Y = 199:GOSUB 25000:NEXT X 
32040 FOR Y=0 TO 199 STEP 6:X=0:GOSUB 
25000:X=319:GOSUB 25000:NEXT Y 

Anche la funzione può essere descritta più rapidamente con un oppor¬ 
tuno STEP, per esempio: 

1060 FOR X=0 TO 319 STEP 3 

Con lo stesso criterio si può aumentare la velocità di tracciamento degli 
assi (righe 10050 e 11050) e degli asintoti (riga 13060). 

Prima di passare alla rassegna di alcuni esempi riteniamo utile segnalare 
una piccola anomalia caratteristica del sistema operativo. 

Con il comando B=INT(A) viene calcolata la parte intera nel numero 
memorizzato nella variabile A. Però, nel caso tale numero sia positivo, 
per esempio A=3.5, si avrà B=3; invece nel caso in cui tale numero sia 
negativo, per esempio A=-3.5, si avrà B=-4. 

Infatti considerando i numeri reali in corrispondenza di una retta orien¬ 
tata verso destra il sistema operativo esegue in entrambi i casi l’appros¬ 
simazione verso sinistra. 


-5 -4 -3 -2 


INT (-3.5) 


INT (3.5) 


Fig. 2.1 


ELENCO DI FUNZIONI DA VISUALIZZARE 


Presentiamo, qui di seguito, una breve raccolta di funzioni di vario tipo 
che il lettore potrà sperimentare scegliendo di volta in volta gli oppor¬ 
tuni intervalli. 
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Funzioni trigonometriche 

(p sta per pi greco) 

y=sin(p/x) 

y=sin(x)/x 

y=tan(x)/(2*sin(x) -1 ) 

y=x*sin(x)/(1 -cos(x)) 

y=(1 -cos(x))/(3*sin(x)*sin(x)) 

y=2*cos(x) - cos(2*x) 

y=(x 12-sin(x))/(cos(x)-1) 

y=(sinx)) 13+(cos(x)) 13 

y=(sin(x/2)) 12 

y=2*tan(x)/(1 +(tan(x)) f 2) 

y=(1 +(tan(x)) 12)/(1 -(tan(x)) f 2) 

y=sin(x)/(1 +cos(x)) 

y=cos(2*x)+3*cos(x) 

y- (sin(x) -2)/(3*sin(x) -1 ) 

y=3*cos(x)/(2*cos(x) -1 ) 

y=(sin(x)-cos(x))/sin(x) 

y=2*(sin(x)) f 2-3*cos(x) 

y=sin(1/x) 

y=(1+cos(x))/sin(x) 


Funzioni logaritmiche 

y=log(1+x)/x 

y=log(x)/x 

y~(x 12-2)/(x*log(x)) 

y=log(abs((x+1)/(x)) 

y=(1 +log(x))/x 

y=2*x-x*log(abs(x)) 

y=log(abs((x+1)/(x-1))) 

y=x*log(abs(x)) 

y=log(abs(x f 2-x-2)) 

y=log(abs((x+3)/(1 -x))) 

y=x+log(abs(2-exp(x))) 

y=log(abs(sin(x))) 



STAMPA DI UNA FUNZIONE CON IL C64 


27 


Funzioni algebriche 

y=(x f 2-5*x+6)/(6-5*x) 

y=(6-5*x)/(x 12-5*x+6) 

y=(x|2+x+1)/(x-1) 

y=(-2*x f 3+3*x f 2+1)/(x f y+2*x-1) 

y=(-3*x+2)/(x f 2-3*x+2) 

y=(4*x-5)/(x 12-1 ) 

y=(x f 2-2*w+1 )/(x t 2+5*x+4) 

y-(x f 3+1)/(x+2) 

y=27/(x 12)-8/(x+1 ) 12 

y=2*x 13/(2*x-1) | 2 

y=xf 2+1/(1-x|2) 

y=(x f 3-3*x+1)/(x f 2-x) 

y=x f 2+x+1)/(2*x 12+3*x+1) 

y=x+3/x-1/x f 2 

y=(x f 3-3*x 12)/(4*x 12-8) 

y=2*x f 2-abs(5*x+2) 

y=(x+3)/abs(x f 2+x-2) 

y=abs(x-1)/(x-1) 

y=abs((2*x 12-3*x)/(x-2)) 

y=abs(-2*x 13+3*x f 2+1 )/(x f 3+2*x-1 ) 

y=x f 3-abs(x) 

y=2*x/(1 +abs(x)) 

y=abs(15*xt5-xt3) 


Funzioni varie 

y=abs(x)*sqr(abs(x))/x 
y=sqr(x f 4-5*x) 
y=1/(x+sqr(1 +x \ 2)) 
y=sqr(x 12-4*x+3)/(x+1) 
y=x+sqr(abs(x f 2-2*x)) 
y=x-sqr(abs(x 12—2*x)) 
y=sqr(abs(x))+1/x 
y=1/x+sqr(abs(1 -x) 12) 
y=x*sqr(abs(x/(x-2))) 
y=(5*x-)/(3*x | (1/3)) 
y=(x+sqr(x))/x 
y=(x+1)*sin(x)/(x 12+2) 
y=x 12*sin(1/x) 
y=x*sin(1/x) 



28 


STAMPA DI UNA FUNZIONE CON IL C64 


y=(2*x f 2+3*x-2)*tan(p*x) 

y=(x f 2-sin(x))/(cos(x)-1) 

y=(x+sin(x))/x 

y=x 12-(sin(x)) f 3+cos(x) 

y=x+2*sin(x) 

y=x*cos(x) 

y=1/(x*cos(x)) 

y=exp(x)*cos(x) 

y=x*cos(x)/2 - si n(x)/2 

y=exp(-x)*sin(x) 

y=*x-(tan(x)) f 2 

y=(x-exp(log(x)))/(x-exp(1 )) 

y=x-sin(x) 

y=x*cos(x)-Sin(x)+p/2 
y=exp( - x)*(sin(x) - 5*cos(x)) 
y=log(abs(1 +x))/(2*sin(x)) 
y=abs(x)*(1 +(sin(1/x)) f 2) 
y=x-3*sin(x/2) 


ESEMPI DI UTILIZZAZIONE 


Dalle figure seguenti si può rilevare l’ottimo grado di definizione con cui 
vengono tracciate le funzioni. Esse rappresentano solo un esempio del 
valido ausilio didattico che il programma consente di ottenere per 
studenti di ogni categoria. 

Negli ultimi due grafici si può constatare quanto già accennato in 
precedenza a proposito delPanomalia legata al comando INT: infatti i 
due grafici dovrebbero essere simmetrici rispetto all’asse y. 
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GRAFICO DELLA FUNZIONE 

V = i :-'*X !2-4*X+5 > / < X T2+1 ) 

ESTREMI DI VISIJAL1 ZZrtZIONE 

x (sinistrh)=-iù 

X < DESTRA>= a 

V (INFERIORE>=-.b 

V (SUPERI ÓRE >= 8 

Fig. 2.2 



GRAFICO DELLA FUNZIONE 

V = CXT2+1 >/(X-XT2> 

ESTREMI DI VISUALIZZAZIONE 

X (SINISTRO ‘=-8 
X (DESTRA)» 4 

V ( INFERIORE "=-4 

V (SUPERIORE)» fc 


Fig. 2.3 
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GRAFICO DELLA FUNZIONE : 

V « X/0+XT4) 

ESTREMI DI VISUALIZZAZIONE 

X (SINISTRA >—5 
X (DESTRA • = A 

V (INFERIORE—? 

V (SUPERIORE)* 2 


Fig. 2.4 



GRAFICO DELLA FUNZIONE - 


V » (<X-1)12)*((X-4)T2> 

ESTREMI DI VISUALIZZAZIONE : 

X (SINISTRA)—.5 
X (DESTRA)* 3 

V (INFERIORE)—.5 

V (SUPERIORE)* 10 


Fig. 2.5 
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GRAFICO DELLA FUNZIONE 

V * 3*S1N(X)-S0R(3>*C0S(X) 

ESTREMI DI VISUÀLIZZA2I0NE : 

X (SINISTRA)—5 
X (DESTRA)* 15 

V (INFERIORE) —IO 

V (SUPERIORE)» IO 


Fig. 2.6 



GRAFICO DELLA FUNZIONE : 


V = 3*$IN(2»W*X)*2»SIN(4*«*X> 


ESTREMI DI VISUALIZZAZIONE • 

X (SINISTRA)*-.5 
X (DESTRA)* ? 

V (INFERIORE)—6 

Y (SUPERIORE)* 6 


Fig. 2.7 
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GRAFICO DELLA FUNZIONE : 

V * SIN'2*;O+2*C0S(X> 
ESTREMI DI VISUALIZZAZIONE 

X (SINISTRA)»- 5 

x «destra:-» le 

V (INFERIORE)»-S 

V (SUPERIORE '» 

Fig. 2.8 



GRAFICO DELLA FUNZIONE : 

V » SINOO/X 

ESTREMI DI VISUALIZZAZIONE : 

x (Sihistra;'=-2ò 
X (DESTRA)» 20 

V (INFERIORE)*-! 

V (SUPERIORE)» 2 


Fig. 2.9 




STAMPA DI UNA FUNZIONE CON IL C64 


33 



GRAFICO DELLA FUNZIONE : 


V • (Xt4*2»Xt3)/(Xt3-l> 

ESTREMI DI VISUALIZZAZIONE : 

X (SINISTRA)—3 
X (DESTRA)- 4 

V (INFERIORE)—4 

V (SUPERIORE)- 8 


Fig. 2.10 



GRAFICO DELLA FUNZIONE : 

V - XT2*SIN(X)/(XT2-1> 

ESTREMI DI VISUALIZZAZIONE : 

X (SINISTRA) —18 
X (DESTRA)- 10 

V (INFERIORE)—10 

V (SUPERIORE)- 10 


Fig. 2.11 
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GRAFICO DELLA FUNZIONE : 

V * X12-ABS < 54X+6) 

ESTREMI DI VISUALIZZAZIONE : 

X (SIHISTRM)=-é> 

X ( DESTRA)» Ò 

V (INFERIORE>*-15 

V (SUPERIORE)* 5 


Fig. 2.12 



GRAFICO DELLA FUNZIONE ' 

V » ABS(X)+ABS(XT2-X) 

ESTREMI DI VISUALIZZAZIONE : 

X (SINISTRA)—4 
X (DESTRA)* 4 

V (INFERIORE)—.5 

V (SUPERIORE)- 20 


Fig. 2.13 
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GRAFICO DELLA FUNZIONE : 

V - SQR(ABS((-l+XT3)/(X+2)>> 

ESTREMI DI VISUALIZZAZIONE : 

X ‘SINISTRA)*-10 
X «DESTRA)* 10 

V (INFERIORE)- 0 

V (SUPERIORE)* 15 


Fig. 2.14 



GRAFICO DELLA FUNZIONE : 


V - ABS(X*LOG(ABS(X>)> 

ESTREMI DI VISUALIZZAZIONE : 

X (SINISTRA)—? 

X (DESTRA)- ? 

V (INFERIORE)—.5 

V (SUPERIORE)* 2 


Fig. 2.15 
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GRAFICO DELLA FUNZIONE ■ 

V * X/INTQO 

ESTREMI DI VISUALIZZAZIONE : 

X (SINISTRA)- 1 
X (DESTRA)- 6 

V (INFERIORE) —.5 

V (SUPERIORE)» 2 


Fig. 2.16 



ORAFICO DELLA FUNZIONE 


V » EXP((X-1 )/(X-H ) ) 

ESTREMI DI VISUALIZZAZIONE : 

X (SINISTRA)—5 
X (DESTRA)» 5 

V (INFERIORE) —.3 

V (SUPERIORE)- 8 


Fig. 2.17 
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GRAFICO DELLA FUNZIONE : 

V - COS(EXP(X>) 

ESTREMI DI VISUALIZZAZIONE : 

X (SINISTRA)*-1 
X (DESTRA)- 3.3 

V (INFERIORE)—3 

V (SUPERIORE)- 3 


Fig. 2.18 



GRAFICO DELLA FUNZIONE 


V - SIN(3*X)+X/2 


ESTREMI DI VISUALIZZAZIONE 

X (SINISTRA)—*5 
X < DESTRA- 1 - 15 

V < INFERIORE) —IO 

V (SUPERIORE>= 1* 


Fig. 2.19 
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GRAFICO DELLA FUNZIONE : 

V » INT(ABS(X))-SOR(l-(X-lNT(X))12) 

ESTREMI DI VISUALIZZAZIONE : 

X (SINISTRA)—5 
X < DESTRA)* 5 

V (INFERIORE)—2 

V (SUPERIORE)* 6 


Fig. 2.20 



GRAFICO DELLA FUNZIONE ■ 

V - INT(ABS(X))-KX-INT(X)>T2 

ESTREMI DI VISUALIZZAZIONE 

X (SINISTRA)—5 
X (DESTRA)» 5 

V (INFERIORE)—2 

V (SUPERIORE)» 6 


Fig. 2.21 





_CAPITOLO TRE 

Stampa di funzioni polidrome 
o sovrapposte con il C64 


DESCRIZIONE DEL LISTATO 


Con i programmi finora realizzati non è possibile visualizzare funzioni 
che abbiano più di una ordinata per ciascuna ascissa (per esempio 
circonferenze, ellissi, iperboli e curve di ordine superiore) oppure più 
funzioni sovrapposte sullo stesso piano cartesiano. 

Pur mantenendo inalterate la struttura di elaborazione, calcolo, trac¬ 
ciamento e stampa, con il seguente programma potremo visualizzare 
quasi ogni tipo di funzione o più funzioni sovrapposte. L’introduzione 
delle funzioni non è più realizzata automaticamente, ma deve essere 
eseguita manualmente inserendo le funzioni stesse nel listato a partire 
dalla riga 440 fino alla riga 570. 

Dopo ogni funzione è necessario aggiungere il comando 
:GOSUB 930 

per il tracciamento. Infatti nell’intervallo 430-590 si ha un ciclo FOR- 
...NEXT con K che va da 0 a 319, durante il quale per ogni ascissa viene 
calcolata e tracciata l’ordinata di ciascuna delle funzioni proposte. 

La riga 580 apparentemente inutile rappresenta infatti un comodo 
espediente, come vedremo fra poco. Per battere correttamente tale 
linea si seguano alla lettera le istruzioni seguenti: 

- battere 580 REM " 

- con il tasto DELETE cancellare le " e ribatterle 
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- battere otto volte il tasto INSERT (mantenendo premuto lo SHIFT) 

- battere altre otto volte il tasto DELETE (senza lo SHIFT) 

- battere RUN 200 

Se ora provate a chiedere LIST 580 noterete sullo schermo la presenza 
della scritta RUN 200 senza il relativo numero di riga: basterà andarci 
sopra con il cursore e premere RETURN per dare il via al programma. 
Possiamo ora capire che grazie a tale riga, dopo aver introdotto le 
funzioni nel listato, possiamo far ripartire il programma senza che sia 
necessario ribattere il comando RUN ed evitando una inutile rilettura 
dei DATA. 

Prima di visualizzare le funzioni è necessario però fissare la scelta degli 
intervalli. 

La stampa si ottiene con il tasto funzione F7 battuto alla fine della 
visualizzazione. 

Ben poco rimane da dire sul resto del programma, molto simile ai 
precedenti. 


100 REM **##***##********************** 


110 REM * COPVRIGHT IN BASIC * 

120 REM * VIA SEZ2E 22 - LATINA * 

130 REM * TEL. 0773 7487631 * 

140 REM *-* 

150 REM * STAMPA FUNZIONI POLIVOGHE * 

160 REM * 0 SOVRAPPOSTE * 


170 REM «Di#*#************##»!********#** 

180 PRINT“CCLEARHGIALL03" : P0KE 53280,14:P0KE 53281,6 

190 GOSUB 1370 

200 PRINT"CCLEARICGIALL0I" 

210 PRINTTAB<5)"t2 DOWNIFt - F7 VISUALIZ. & STAMPA" 

220 PRINTTAB<5)"F3 - CAMBIO DELLE FUNZIONI" 

230 PRINTTAB<5>”F5 - SCELTA DEGLI INTERVALLI" 

240 GOSUB 850 

250 IF QGC133 OR QQ>135 THEN 240 

260 0N QG-132 GOTO 270,410,610 

270 REM - 

280 REM VISUALIZZAZIONE 

290 REM - 

300 SVS 49152 

310 DX=»<X2-X1 >7320 : DV=<Y2-Y1 >7200 

320 GOSUB 980:GOSUB 1090 REM ASSI COORDINATI 

330 GOSUB 430 

340 GOSUB 850 

350 IF QQ=136 THEN GOSUB 780 

360 GOSUB 1280 
370 GOTO 200 

380 REM - 

390 REM CAMBIO FUNZIONI 

400 REM - 

410 PRINT"CH0MEK6 D0WNI" 

420 LIST 450-580 

430 F0R K=0 T0 319 

440 X=X1+K*DX : GOSUB 1140 

450 V=SIN<X72>*C0SCX> : GOSUB 930 

460 Y*172*SIN<X>#C0S<2#X>:GOSUB 930 

470 Y=173#SIN(1.5#X)*C0S(3#X> : GOSUB 930 
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480 V*SIN<X/2)+lX2#SIN(X)+lX3#SIN(1.5*X) : G0SUB 930 
490 : 

500 : 

510 ^ 

520 ^ 

530 ^ 

540 : 

550 ■ 

560 : 

570 • 

580 REM "iMMMIIdJN 200 
590 NEXT K 
600 RETURN 

610 REM - 

620 REM CAMBIO INTERVALLI 


Doti r.cn - _ 

640 PRINT"CCLEAR3SCEGLI GLI INTERVALLI DI VISUALIZZAZIONE 

650 INPUT "12 DOWN3X (SINISTRA)*";XI 

660 INPUT "X (DESTRA)*";X2 

670 IP X1>=X2 THEN 640 

680 INPUT "CDOWN3V (BASSA)*";VI 

690 INPUT "V (ALTA)*";V2 

700 IF VI >=V2 THEN 680 

710 PRINT"CDOWN3CI SONO ASINTOTI VERTICALI ?" i 
7£0 orisiIB 850 PRINTQi IF Q#*"N" THEN NA=-l-’GOTO 770 
730 FRI NT" C DOWN 3 QUANTI GOSUB 850 : F'RINTQ 

740 NA=G FOR K=1 TO HA 

750 PRINT"CDOWN3ASCISSA"K"* ASINTOTO = "; : INPUT AS(K) 

760 NEXT K 
776 GOTO 200 


780 REM - 

790 REM STAMPA FUNZIONI 

800 REM - 

810 REM - TRACCIAM. CORNICE 


820 FOR K=0 TO 319 STEP 3:V=0 GOSUB 1210:V=199 :GOSUB 12 

10 NEXT K 

830 FOR Y=0 TO 199 STEP 3 : K=8 GOSUB 1210 ■' K=319 : GOSUB 12 

10 NEXT V 

840 SVS 52224 RETURN 

850 REM - 

860 REM LOOP DI ATTESA 

870 REM - 

880 OET 0$ '■ IF Qt="" THEN 880 

A?fl QQ=ASC<0#> : Q=VAL( OS"> RETURN 

900 REM - 

910 REM ELABORAZIONE COORDINATE 

920 REM - 

930 IF YOY1 THEN 970 
940 IF V>=V2 THEN 970 
950 V=200-<V-V1>/DV 
960 GOSUB 1210 
970 RETURN 

980 REM - 

990 REM ASSE X 

1000 REM - 

1010 IF SGN(VI)* SGN(V2) THEN RETURN 

1020 V*INT(V2/BV+.5) 

1030 FOR K=0 TO 319 STEP 6 
1040 GOSUB 1210 
1050 NEXT K RETURN 

1060 REM - 

1070 REM ASSE V 

1080 REM - 

1090 IF SGN(Xl)* SGN(X2) THEN RETURN 
1100 K=-INT(X1/DX+,5) 
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1110 FOR V=0 TO 199 STEP 6 

1120 GOSUB 1210 
1130 NEXT V : RETURN 

1140 REM - 

1150 REM ASINTOTI 

neo REM - 

1165 IF NA=-1 THEN RETURN 

1170 FOR J=*l TO NfUIF X=AS<J) THEN 1190 

1180 NEXT J : RETURN 

1190 FOR T=0 TO 199 STEP 6:V=T:G0SUB 1210:NEXT T 

1200 X=X1+<K+1)#DX:RETURN 

1210 REM - 

1220 REM PLOTTAMENTO 

1230 REM - 

1240 IF K<0 OR K>319 OR Y<0 OR V>199 THEN 1270 


1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 


BV=8192+320*INT<Y/3)+8*1NT<K/8)+<V AND 7) 

POKE BV,PEEK<BY> OR 2t<7-<« AND 7)> 

RETURN 

REM - 

REM DISATTIVAZIONE COMANDO PLOT 

REM - 

POKE 53272,PEEK<53272) AND 247 
POKE 53265,PEEK<53265) AND 223 
PRINT"[OLEARI"'RETURN 

REM - 

REM ROUTINE LETTURA DATA 

REM - 

PRINTTAB<10)"C2 DOWN]ATTENDERE PREGO _" 

FOR K- 49152 TO 49222 - 

READ NPOKE K,N : ND=ND+N : NEXT K 

POKE 2, 2 -REM PER STAMPA RISTRETTA SOSTITUIRE CON POKE 2, 
1 

FOR K=52224 TO 52677 

READ N : POKE K,N : ND=ND+N : NEXT K 

IF NDO70923 THEN PRINT 11 CDOWN)ERRORE NEI DATA" :END 
RETURN 

DATA 173,24,203,9,8,141,24,208 
DATA 173,17,208,9,32,141,17,208 
DATA 169,0,133,251,169,32,133,252 
DATA 160,0,169,0,145,251,200, 192 
DATA 0,208,249,230,252,169,64,197 
DATA 252,208,239,169,0,133,251,169 
DATA 4,133,252,160,0,169,3,145 
DATA 251,200,192,0,208,249,230,252 
DATA 169,8,197,252,208,239,96 

REM - 

REM ROUTINE PER HARD COPV 

REM - 

DATA 173,0,221,41,3,73 
DATA 3,160,6,10,136,208 
DATA 252,141,132,3,173,24 
DATA 208,41,8,240,9,24 
DATA 169,32,109,132,3,141 
DATA 132,3,169,0,32,189 
DATA 255,169,4,170,160,255 
DATA 32,186,255,32,192,255 
DATA 162,4,32,201,255,176 
DATA 3,76,61,204,76,32 
DATA 205,169,8,32,210,255 
DATA 169,13,32,210,255,162 
DATA 0,169,1,141,198,205 
DATA 169,0,141,199,205,169 
DATA 0,141,200,205,169,199 
DATA 141,201,205,32,225,255 
DATA 208,3,76,32,205,138 
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1740 DATA 
1750 DATA 
1760 DATA 
1770 DATA 
1780 DATA 
1790 DATA 
1800 DATA 
1810 DATA 
1820 DATA 
1830 DATA 
1840 DATA 
1850 DATA 
1860 DATA 
1870 DATA 
1880 DATA 
1890 DATA 
1900 DATA 
1910 DATA 
1920 DATA 
1930 DATA 
1940 DATA 
1950 DATA 
1960 DATA 
1970 DATA 
1980 DATA 
1990 DATA 
2000 DATA 
2810 DATA 
2020 DATA 
2030 DATA 
2040 DATA 
2050 DATA 
2060 DATA 
2070 DATA 
2080 DATA 
2890 DATA 
2100 DATA 
2110 DATA 
2120 DATA 
2130 DATA 
2140 DATA 
2158 DATA 
2160 DATA 
2170 DATA 
2180 DATA 
2190 DATA 
2200 DATA 
2210 DATA 
2220 DATA 
2230 DATA 
2240 DATA 
2250 DATA 
2260 DATA 
2270 DATA 
2280 DATA 
2290 DATA 
2300 DATA 
2310 DATA 
2320 DATA 


72.152.72.32.41.205 

104.168.104.170.173.205 
205,45,203,205,240,12 

173.202.205.13.198.205 

141.202.205.76.143.204 

173.198.205.73.255.45 
202,205,141,202,205,14 
193,205,173,198,205,201 

128.240.20.24.173.199 

205.105.1.141.199.205 
173,200,205,105,0,141 

200.205.76.93.204.173 

202.205.9.128.224.45 

144.10.173.202.205.41 

31.9.128.141.202.205 
142,207,205,166,2,142 
206,205,168,32,210,255 
152,206,206,285,288,246 

174.207.205.169.1.141 
198,205,169,0,141,202 
205,56,173,199,205,233 

6.141.199.205.173.200 
205,233,0,141,200,205- 

206.201.205.173.201.205 
201,255,240,3,76,93 
204,224,45,176,31,24 

173.199.205.105.7.141 
199,205,173,200,205,105 
0,141,200,205,232,169 
199,141,201,205,169,13 

32.210.255.76.93.204 
169,13,32,210,255,32 
231,255,96•173,201,205 

41.7.141.204.205.173 
201,205,74,74,74,168 
185,146,205,133,251,185 
172,205,138,252,24,165 
251,109,204,205,133,251 
165,252,105,0,133,252 
24,173,132,3,101,252 

133.252.173.199.205.41 
7,73,7,168,200,169 
0,56,42,136,208,252 

141.203.205.24.173.200 

205.101.252.133.252.173 
199,205,41,248,168,138 

72.120.162.52.134.1 

177.251.162.55.134.1 

88.168.104.170.152.45 
203,285,141,205,205,96 
0,64,128,192,0,64 

128,192,0,64,128,192 

0,64,128,192,0,64 

128,192,0,64,128,192 

0,64,0,1,2,3 

5,6,7,8,10,11 

12,13,15,16,17,18 

20,21,22,23,25,26 

27,28,30,31 
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ESEMPI DI UTILIZZAZIONE 


Nella figura seguente si può notare il tracciamento di una parabola con 
asse obliquo sovrapposto a una funzione trascendente. La descrizione 
della parabola 


x 2 - 2xy + y 2 - x - 2 = 0 


si realizza spezzandola in due rami corrispondenti alle due funzioni che 
si ottengono risolvendo l’equazione rispetto alla variabile y: 

y 2 - 2xy + (x 2 - x - 2) = 0 



ESTREMI DI VISUALIZZAZIONE : 

X (SINISTRA)»-2.5 
X (DESTRA)» 4 

V (INFEPI0RE>»-2.3 

V (SUPERIORE)» 5 


450 V»EXP < X >/X T2 : GOSUB 930 
460 IF X<-2 THEN 590 
470 V»X-SQR(X+2):GOSUB 930 
480 V-X+SQRCX+2):GOSUB 930 
490 ^ 

500 : 

510 : 

520 
530 : 

540 : 

550 : 

560 : 

570 : 


Fig. 3.1 
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Applicando la formula risolutiva ridotta si ha 
y - x ± Jx + 2 

Con il segno positivo si avrà il ramo superiore mentre con il segno 
negativo quello inferiore. 

Il tracciamento di questi due rami si verifica alle righe 470 - 480 mentre, 
per evitare che valori della x minori di -2 rendano negativo il radicando 
provocando un segnale di errore, nella riga precedente (460) viene 
imposto un salto condizionato che impedisce il passaggio alle linee 
sottostanti. 

Nella prima riga, la 450, viene definita invece una funzione trascenden¬ 
te che sarà invece tracciata integralmente in quanto si trova a monte del 
salto condizionato. 

Si ricordi di notificare al calcolatore resistenza dell’asintoto verticale 
per x = 0. 



ESTREMI DI VISUALIZZAZIONE : 

X (SINISTRA)—3 
X (DESTRA)- 3 

V (INFERIORE)—10 

V (SUPERIORE)- 10 


450 V-(X-2)*EXP(X+1> QOSUB 930 
460 V-EXP(1-COS(2*X)) : GOSUB 930 
470 : 

488 : 

490 ■ 

300 : 

310 ^ 

320 : 

330 : 

340 : 

330 ■ 

360 : 

370 : 


Fig. 3.2 







46 


STAMPA DI FUNZIONI POLIDROME O SOVRAPPOSTE CON IL C64 


Nessuna precauzione particolare è invece necessaria per ottenere la Fig. 
3.2. La Fig. 3.3 mostra invece un simpatico fiocco ottenuto al solito 
suddividendo la funzione 

y 2 = x 3 + x 2 

nei due rami 

y = ± Jx^Tx 2 

con la solita condizione per evitare il calcolo nei casi in cui il radicando 
diverrebbe negativo. 



ESTREMI DI VISURLIZZfiZIONE • 

X (SINISTRA)—2 
X (DESTRA)- 2 

V (INFERIORE)—2 

V (SUPERIORE)" 2 


430 V«X:GOSUB 930 

460 V—X : GOSUB 930 

470 IF X<-1 THEN 390 

480 V-SQR(XT3+X12 ):GOSUB 930 

490 V—SQR(XT3+Xt2r00SUB 930 

300 : 

310 : 

320 : 

330 : 

340 : 

550 : 

560 : 

370 : 


Fig. 3.3 
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Nelle ultime due figure si ha invece una combinazione di fantasia fra 
coniche e funzioni trigonometriche. 

Si noti che la circonferenza viene riprodotta come un’ellisse perché sui 
due assi si hanno unità di misura differenti; infatti pur essendo identici 
gii estremi di visualizzazione orizzontale e verticale, i pixel disponibili 
per l’asse x sono 320 mentri quelli disponibili per l’asse y sono soltanto 
200 . 



V 

/ / 



.. J: .X 

/ 



ESTREMI DI VISUALIZZAZIONE ■' 

X (SINISTRA}*-!.3 
X (DESTRA)» 1.3 

V (INFERIORE)=-1.3 

V (SUPERIORE)* 1.5 


430 IF X<-1 OR X>1 THEN 590 

460 V*S0R(1-X t2) : GQSUB 930 

470 V=-SQR(1-Xt2) OOSUB 930 

480 V»SIN(2*fr*X)*SQR(l-Xt2) :G0SUB 930 

490 

300 : 

510 : 

320 : 

530 : 

540 : 

330 : 

560 : 

370 : 


Fig. 3.4 
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ESTREMI DI VISUALIZZAZIONE : 

X (SINISTRA)«-5 
X <DESTRA)® 3 

V (INFERIORE)—3 

V (SUPERIORE)® 5 


450 

460 

470 

480 

490 

500 

510 

320 

530 

540 

550 

560 

570 


V®-X : GOSUB 930 

V=X : GOSUB 930 

IF X>-1 AND X<1 THEN NEXT 

V=SQP(X T2-1):GOSUB 930 

V=-SC'P(Xt2-l) GOSUB 938 

V«C0S(2#X)*S0R(XT2-1) GOSUB 930 


Fig. 3.5 


APPLICAZIONI NELLA SCOMPOSIZIONE DI FOURIER 


Ogni funzione (anche non periodica) può essere scomposta in un solo 
modo in una somma di opportune oscillazioni armoniche 

f(x ) = all + a x x cos(x) + bi x *sin(x) + a 2 x cos(2x) + b 2 x sin(2x) + 

dove i coefficienti sono calcolabili con le formule 

i C 2k 

a k = J f{x) X cos (kx) dx 
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1 C 2n 

b k = — I f(x) x sin ( kx) dx 

n J o 

Tale principio è noto come principio di Fourier. 

Questo teorema è di importanza didattica fondamentale e il programma 
che stiamo trattando si presta egregiamente alla visualizzazione delle 
armoniche e della funzione /( x) risultante dalla loro composizione. 
Può risultare utile ricordare che: 

- se il valore medio di f(x) è nullo allora il termine all è assente; 

- se/(jc) è una curva alternata si hanno solo le armoniche dispari perché 
solo per esse si può ottenere: 

f(t + T/2) = -f(t) 

cioè la curva assume ordinata opposta dopo un semiperiodo; 

- se f(x) è simmetrica rispetto all’asse y allora i termini in seno sono 
nulli perché per tali termini si ha 

sin(A) - -sin(-A) 

e pertanto per essi non può mai avvenire che 
/(')=/(-') 

- s ef(x) è simmetrica rispetto all’origine deve invece essere verificata la 
relazione 

m = -f(-t) 

e ciò può essere verificato solo per i termini in seno in quanto 
cos(A:) = cos(—A:) 

Passiamo all’analisi di alcune funzioni caratteristiche. 


Onda triangolare 

È composta da infiniti termini dispari in seno di ampiezza decrescente e 
di segno alternato. 

y = sin(;t)—1/9 x sin(3^f) + l/25xsin(5jt)—... 
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Per la sua visualizzazione è sufficiente limitarsi alle prime tre armo¬ 
niche. 

L’intervallo di rappresentazione è 

x (sinistra) = 0 x (destra) = 6.28 

y (inferiore) = -2 y (superiore) = 2 



Fig. 3.6 Onda triangolare. 


Onda quadra 

y = sin(x) + 1/3 x sin(3x) + 1/5 x sin(5x) + ... 
Estremi di visualizzazione identici ai precedenti. 
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Onda a dente di sega 

y = jt/2 — sin (x) - 1/2 x sin(2x) -1/3 x sin(3x) - ... 

Intervallo di visualizzazione 

x (sinistra) = 0 x (destra) = 12.56 

y (inferiore) = 1.1 y (superiore) = 3.2 
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Fig. 3.8 Onda a dente di sega. 


FUNZIONI RECIPROCHE 


Una interessante applicazione che proponiamo è costituita dalla possi¬ 
bilità di visualizzare su uno stesso piano cartesiano una generica fun¬ 
zione 

y = m 

e la sua corrispondente reciproca 

y= W*) 

Queste due funzioni sovrapposte possono mostrare interessanti rela¬ 
zioni in quanto godono delle seguenti proprietà: 
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- tracciando le rette orizzontali y=le y=-l la f(x) e la sua reciproca 
andranno a intersecarle sempre in punti coincidenti; 

- ogni intersezione della f(x ) con le ascisse corrisponde a un asintoto 
verticale per la reciproca e viceversa; 

- ogni asintoto orizzontale viene trasformato in un altro asintoto oriz¬ 
zontale con ordinata reciproca; 

- ogni punto di massimo si trasforma in un punto di minimo e viceversa 
(le ordinate saranno ovviamente reciproche). 



ESTREMI PI VISUALIZZAZIONE : 

X (SINlSTRA)»-3 
X (DESTRA)» 6 

V (INFERI0RE>»~6 

V (SUPERIORE)» 8 


438 V»-l OOSUB 938 
468 V»1 : OOSUB 938 
478 V»XT2-5*X+6 : OOSUB 930 
480 V»l/(XT2-3*X+6) OOSUB 930 

490 : 

500 : 

310 : 

320 : 

330 : 

340 : 

330 : 

360 : 

570 : 


Fig. 3.9 
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ESTREMI DI VISUALIZZAZIONE 

X (SINISTRA)—3 
X (DESTRA)* 10 

V (INFERIORE)—5 

V (SUPERIORE)* 13 


450 V—PGOSUB 930 

460 V*1-OOSUB 930 

470 V»(X-2)X(X T2+X-2) : OOSUB 930 

480 V*(XT2+X-2)/(X-2> OOSUB 930 

490 •• 

500 : 

310 : 

320 : 

330 : 

340 : 

330 : 

560 : 

570 : 


Fig. 3.10 



Fig. 3.11 
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ESTREMI DI VISUALIZ2AZIONE : 

X (SINISTRA)«-3.3 
X <DESTRA)* 3.5 

V (INFERIORE)=-4 

V (SUPERIORE)- 6 


450 V—1 : OOSUB 330 
460 V-l : OOSUB 930 
470 V=XT4-5*XT2+4 OOSUB 930 
480 V»l/(Xt4-5*XT2+4) OOSUB 930 
490 : 

500 : 

510 : 

520 : 

530 ^ 

540 : 

550 : 

560 : 

570 ■■ 


Fig. 3.11 (segue). 


SOVRAPPOSIZIONE FRA DUE ONDE: BATTIMENTI 


Un ulteriore aspetto, interessante dal punto di vista didattico, è la 
possibilità di mostrare come due segnali sinusoidali possano combinarsi 
generando un segnale complesso avente ventri e nodi disposti periodi¬ 
camente. 

Prendiamo in considerazione due funzioni sinusoidali aventi per sem¬ 
plicità la stessa ampiezza a 

y x - a sin (vi x) 
y 2 = a sin (v 2 x) 

v 2 e v 2 rappresentano le pulsazioni delle due onde (cioè la loro frequen¬ 
za moltiplicata per 2 it). Sommando le due funzioni si ottiene la nuova 
funzione 

y = y x + y 2 = a (sin (vj x) + sin (v 2 x)) 

Applicando la formula di prostaferesi 

sin (p) + sin ( q) = 2 sin ((p + q)l2 ) cos ((p - q)/2) 

si ha 

y = 2 a cos (x (v t - v 2 )/2) sin (x (v t + v 2 )/2) 
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Questa funzione è formata da un fattore in coseno moltiplicato per un 
secondo fattore in seno; essa è ancora una funzione periodica, ma con 
ampiezza massima dipendente dal tempo e non costante. 



ESTREMI DI VISUALIZZAZIONE : 

X (SINISTRA) —1 
X <DESTRA)* 10 

V < INFERIORE)—5 

V (SUPERIORE)* 5 


450 V«C0S((5-4)*X)*SIN((5+4)*X) GOSUB 930 
460 • 

470 : 

480 ^ 

490 : 

500 : 

510 : 

520 : 

530 ^ 

540 : 

550 : 

560 : 

570 : 


Fig. 3.12 


Il fattore coseno determina l’ampiezza massima dipendente dal tempo 
mentre il fattore seno genera delle oscillazioni più frequenti la cu 
ampiezza è fissata in modo variabile dal fattore precedente. 

La figura mostra chiaramente il segnale sinusoidale modulato in am 
piezza. 
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FIGURE DI LISSAJOU 


Si abbiano due moti armonici 

x = a sin (vj t) 
y = b sin (v 2 t + s) 

dove v t e v 2 rappresentano le pulsazioni, s lo sfasamento tra i due moti e 
a e b le rispettive ampiezze massime. 

Combiniamoli considerandoli come le coordinate parametriche di un 
punto: in altre parole la figura risultante dalla combinazione avrà ascissa 
e ordinata coincidenti con i due moti armonici. Il parametro rè il tempo. 
Su alcuni testi di fisica sono spesso riportate tali figure, denominate 
figure di Lissajou, e ottenute per la prima volta con il doppio pendolo 
del Righi. Lo svolgimento di questo argomento può utilmente essere 
sfruttato per applicazioni didattiche durante la trattazione dei moti 
armonici. 

Il luogo geometrico descritto dal punto avente le suddette coordinate 
parametriche può ottenersi eliminando il parametro t, ricavandolo dalla 
prima funzione e sostituendolo nella seconda. 

Cominciamo ad analizzare il caso in cui vj = v 2 , ricordando che 


cos (vi t) = sin /(I — (sin (vq t)) 2) 

Eseguendo i calcoli si ricava l’equazione di una ellisse ruotata con centro 
nell’origine: infatti la presenza del termine in xy denuncia la rotazione, 
mentre la mancanza dei termini lineari in x e y evidenzia il fatto che il 
centro di simmetria coincide con l’origine degli assi. 

y x : 

2 u -COS (s ) \ 2 

4+ r -‘m - 1 

a \ sin ( s) ) 

Nel caso in cui a = b si ha una circonferenza. 

È interessante la trattazione della evoluzione della figura quando lo 
sfasamento s varia fra 0 e 180 gradi. 

Per s — 0 si ha y = ( b/a)x , cioè un segmento rettilineo con centro 
nell’origine degli assi e coefficiente angolare bla. 

Per valori compresi fra 0 e 90 gradi si ottiene una ellisse con centro 
nell’origine e asse maggiore inclinato verso l’alto. 

Per s = 90 gradi l’ellisse è orizzontale con equazione canonica 
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Per valori compresi fra 90 e 180 gradi si ottiene di nuovo una ellisse, ma 
ruotata verso il basso. 

Infine per s = 180 gradi si ha un segmento rettilineo di equazione 
y = (—b/a)x. 

Se invece V! e v 2 non sono uguali si ottiene una figura variamente 
intrecciata. 


10 REM ******************************** 

20 REM * TRACCI PIMENTO * 

30 REM * FIGURE DI LISSATOLI * 

40 REM ******************************** 

50 CLR : POKE 53230,1 POKE 53231,6 
60 PRINT"COLEARI":GOSUB 1340 

70 PRINT"[OLEARICDOWN3"TAB<10)"FIGURE DI LISSAJOU" 

80 INPUT "CDOWN3PULSAZIONE ORIZZONTALE",PX 
85 INPUT "IDOWNIPULSAZIONE VERTICALE",PV 
90 INPUT "CD0WN3SFASAM. TRA I DUE MOTI <IN GRADI>";S 
100 SVS 49152 
110 K=0S=S*ff/180 
120 X=160*SIN<K*PX>+160 
130 V=200-100*SIN <K*PV+S > +100 
200 GOSUB 1240 

500 K=K+.05 IF K*PX<2*ff OR K*PV+S<2*ff THEN 120 
510 GOSUB 830 

520 IF 0*="+-" THEN GOSUB 820 
530 GOSUB 1310 

540 PRINT"COLEARIANCORA ?" GOSUB 880 

550 IF C'*="S" THEN 70 

640 PRINT"12 DOWNIOK - FINE LAVORO," : END 

780 REM - 

790 REM STAMPA FIGURA 

800 REM - 


810 

REM 


- TR 

ACCIAI' 

1. CORNICE- 



820 

FOR 

K-0 

TO 

319 

STEP 3 :V=0:GOSUB 

1210:y=199:GOSUB 

12 


10NEXT K 






830 

FOR 

V=0 

TO 

199 

STEP 3 :K=0:GOSUB 

1210 : K=319 : GOSUB 

12 


10NEXT V 


840 SVS 52224:RETURN 

350 REM - 

860 REM LOOP DI ATTESA 

870 REM - 

880 GET Q$ : IF Q*="" THEN 880 
890 QQ=ASC<Q*):Q=VAL<Q$): RETURN 

1210 REM - 

1220 REM PLOTTAMENTO 

1230 REM - 

1240 IF X<0 OR X>319 OR V<0 OR V>199 THEN 
1250 BV*3192+320*INT<V/8>+8*1NT<X/8)+<V AND 7) 
1260 POKE BV,PEEK<BV> OR 2t(7-(X AND 7>> 

1270 RETURN 

1280 REM - 

1290 REM DISATTIVAZIONE COMANDO PLOT 

1300 REM - 


1270 
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1310 POKE 53272, PEEK< 53272) AND 247 
1320 POKE 53265, PEEK< 53265) AND 223 
1330 PRINT" l CLEARI"■RETURN 

1340 REM - 

1350 REM ROUTINE LETTURA DATA 

1360 REM - 

1370 PRINTTRB<10)"12 DOWNIATTENDERE PREGO 

1380 POR K* 49152 TO 49222 

1390 READ N : POKE K,N:ND=ND+N:NEXT K 

1400 POKE 2 ,2 REM PER STAMPA RISTRETTA SOSTITUIRE CON POKE 2, 
1 

1410 FOR K*52224 TO 52677 

1420 READ N : POKE K,N:ND=ND+N:NEXT K 

1430 IF NDO70923 THEN PRINT"CDOWNTERRORE NEI DATA": END 
1440 RETURN 

1450 DATA 173,24,208,9,8,141,24,208 
1460 DATA 173,17,208,9,32,141,17,208 
1470 DATA 169,0,133,251,169,32,133,252 
1480 DATA 160,0,169,0,145,251,200,192 
1490 DATA 0,208,249,230,252,169,64,197 
1500 DATA 252.208,239,169,0,133,251,169 
1510 DATA 4,133,252,160,0,169,3,145 
1520 DATA 251,200,192,0,208,249,230,252 
1530 DATA 169,8,197,252,208,239,96 

1540 REM - 

1550 REM ROUTINE PER HARD COPV 

1560 REM - 

1570 DATA 173,0,221,41,3,73 
1580 DATA 3,160,6,10,136,208 
1590 DATA 252,141,132,3,173,24 
1600 DATA 208,41,8.240,9,24 
1610 DATA 169,32,109,132,3,141 
1620 DATA 132,3,169,0,32,189 
1630 DATA 255,169,4,170,160,255 
1640 DATA 32,186,255,32,192,255 
1650 DATA 162,4,32,201,255,176 
1660 DATA 3,76,61,284,76,32 
1670 DATA 205,169,8,32,210,255 
1680 DATA 169,13,32,218,255,162 
1690 DATA 0,169,1,141,198,205 
1700 DATA 169,0,141,199,285,169 
1710 DATA 0,141,200,285,169.199 
1720 DATA 141,201,205,32,225,255 
1730 DATA 283,3,76,32,205,138 
1740 DATA 72,152,72,32,41,205 
1750 DATA 104,168,104,170,173,285 
1768 DATA 205,45,203,205,240,12 
1770 DATA 173,202,205,13,198,205 
1780 DATA 141,202,205,76,143,204 
1790 DATA 173,198,285,73,255,45 
1800 DATA 202,205,141,202,285,14 
1810 DATA 193,205,173,198,205,201 
1820 DATA 128,240,20,24,173,199 
1330 DATA 205,105,1,141,199,205 
1840 DATA 173,200,205,105,0,141 
1850 DATA 200,205,76,93,204,173 
1868 DATA 282,205,9,128,224,45 
1870 DATA 144,10,173,202,205,41 
1880 DATA 31,9,128,141,202,205 
1890 DATA 142,207,205,166,2,142 
1900 DATA 286,205,168,32,210,255 
1910 DATA 152,206,206,205,288,246 
1920 DATA 174,207,205,169,1,141 
1938 DATA 198,205,169,0,141,202 
1940 DATA 205,56,173,199,205,233 
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1950 DATA 
1960 DATA 
1970 DATA 
1980 DATA 
1990 DATA 
2000 DATA 
2010 DATA 
2020 DATA 
2030 DATA 
20-40 DATA 
2050 DATA 
2060 DATA 
2070 DATA 
2080 DATA 
2090 DATA 
2100 DATA 
2110 DATA 
2120 DATA 
2130 DATA 
2140 DATA 
2150 DATA 
2160 DATA 
2170 DATA 
2180 DATA 
2190 DATA 
2200 DATA 
2210 DATA 
2228 DATA 
2230 DATA 
2240 DATA 
2250 DATA 
2260 DATA 
2270 DATA 
2280 DATA 
2290 DATA 
2300 DATA 
2310 DATA 
2320 DATA 


6.141.199.205.173.200 
205,233,0,141,200,205 

206.201.205.173.201.205 
201,255,240,3,76,93 
204,224,45,176,31,24 
173-199,205,105,7,141 
199,205,173,200,205,105 
0,141,200,205,232,169 
199,141,201,205,169,13 
32,210,255,76,93,204 
169,13,32,210,255,32 

231.255.96.173.201.205 

41.7.141.204.205.173 
201,205,74,74,74,168 
185,146,205,133,251,185 
172,205,133,252,24,165 
251,109,204,205,133,251 
165,252,105,0,133,252 
24,173,132,3,101,252 
133,252,173,199-205,41 
7,73,7,168,200,169 
0,56,42,136,208,252 

141.283.205.24.173.200 

205.101.252.133.252.173 
199,205,41,248,168, 138 
72,120,162,52, 134, 1 
177,251,162,55,134,1 
88,168,104,170,152,45 
203,205,141,205,205,96 
0,64,128,192,8,64 

128,192,0,64,128,192 
0,64,128-192,0,64 
128,192,0,64,128,192 
0,64,0,1.2,3 
5,6,7,8,10,11 
12,13,15,16,17,18 
20,21,22,23,25,26 
27,23,30,31 


Le pulsazioni Vj e v 2 corrispondono nel listalo alle variabili PX e PY, le 
ampiezze a e b corrispondono rispettivamente ai valori 160 e 100 (in 
modo che la figura risultante venga a tutto schermo). 

Le due funzioni sono introdotte nelle righe 120 e 130 e si noti come è 
necessario aggiungere alle coordinate i valori 100 e 160 per trasportare 
l’origine degli assi nel punto centrale dello schermo. Inoltre per le 
ordinate è necessario determinare il complemento a 200 per ovviare al 
fatto che l’asse y dello schermo è orientato verso il basso, mentre nel 
piano cartesiano esso è normalmente orientato verso l’alto. 

Il parametro t nel nostro caso è rappresentato dalla variabile K che, nella 
riga 500, viene successivamente incrementata fino a quando almeno 
una delle due componenti non abbia eseguito un ciclo completo. 
Poiché l’argomento del seno deve essere espresso in radianti, nella riga 
110 lo sfasamento viene trasformato da gradi in radianti. 

Per il resto il programma segue fedelmente i criteri già visti. 
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Fig. 3.13 


Nella prima figura abbiamo posto PX = PY = 1, S = 45 gradi, ottenen¬ 
do così un’ellisse ruotata verso l’alto come previsto teoricamente. 
Nella seconda figura invece è stato posto PX = 1, PY = .5, S = 0. 
Infine nella terza figura invece è stato posto PX = 1, PY = 3, S = 90, ed 
è stato necessario trasformare l’incremento .05 alla riga 500 in un 
incremento più piccolo .01, altrimenti i punti della figura sarebbero 
risultati troppo distanziati fra loro. 





62 


STAMPA DI FUNZIONI POUDROME O SOVRAPPOSTE CON IL C64 




Fig. 3.14 
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Fig. 3.15 









CAPITOLO QUATTRO 

Stampa di funzioni 
con il Simon’s Basic 


DESCRIZIONE DEL LISTATO 


Ci occuperemo ora di un programma che richiede l’impiego del Simon’s 
Basic, un importante strumento che aggiunge molti nuovi comandi 
(oltre 100) al Basic standard del Commodore 64. 

È un programma in linguaggio macchina (facilmente reperibile, ma che 
non possiamo fornire in questo libro perché coperto da Copyright) che 
va caricato e al quale va dato il RUN. 

Da quel momento il Commodore 64 è in grado, con opportuni e semplici 
comandi Basic, di operare in alta risoluzione, di avere una gestione 
facilitata del suono, di trasmettere su stampante il contenuto dello 
schermo, di operare con programmazione strutturata ecc. 

Noi ci limiteremo a usare soltanto alcuni comandi di tracciamento e 
Yhardcopy, che ci eviteranno il ricorso alle routine in DATA e inoltre 
consentiranno una esecuzione dei grafici molto più veloce. 

In più sarà possibile scrivere con caratteri grafici sullo schermo in alta 
risoluzione. 


160 REM ******************************* 
110 REM * COPVRIGHT IN BASIC * 
120 REM * VIA SEZ2E 22 - LATINA * 
130 REM * TEL. 0773 7487631 * 

140 REM *-* 

150 REM * STAMPA FUNZIONI POLIVOGHE * 
160 REM * 0 SOVRAPPOSTE * 
170 REM Ì»!##*!#***###*#**####*##**##**### 




STAMPA DI FUNZIONI CON IL SIMON’S BASIC 


65 


180 PRINT'TJil" ' POKE 53280» 14 : POKE 53281,6 
21A PRINTTAB( 5 > ”KUF 1 - F7 VISUALIZ. fi STAMPA" 
220 PRINTTflB(5)"F3 - CAMBIO DELLE FUNZIONI" 
230 PRTNTTAB(5)"F5 - SCELTA DEGLI INTERVALLI" 
240 GQSUB 850 

250 IF QGC133 OR Q0>135 THEH 240 
260 OH OQ-132 GOTO 270,410,610 

270 REM - 

280 REM VISUALIZZAZIONE 

290 REM - 

380 HIRES 6,3 

310 DX= < X2-X1 ) /320 : DV= < V2-V1 ) /200 

320 OOSUB 980 GOSUB 1030 : REM ASSI COORDINATI 

330 GOSUB 430 

340 GOSUB 850 

350 IF QQ=136 THEH GOSUB 780 
970 HRM GOTO 180 

380 REM - 

390 REM CAMBIO FUNZIONI 

4O0 REM - 

410 PRINT"«W«*m»W" 

420 LIST 450-580 

430 FOR K=0 TO 319 

440 X=X1+K*DX ; GOSUB 1140 

450 y=-l GOSUB 930 

460 V=1 GOSUB 930 

470 Y=(X-2V<XT2+X-2) : GOSUB 930 

480 V= ( X T2+X-2 >/<X-2) : GOSUB 930 

490 • 

500 : 

510 : 

520 ^ 

530 : 

540 •' 

550 - 
560 : 

570 ^ 


580 REM"lta*«*a*F:UN 
590 HEXT K 
600 RETURN 

610 REM - 

620 REM CAMBIO INTERVALLI 


DOv r.c.1 i- 

640 PRINT”3SCEGLI GLI INTERVALLI DI VISUALIZZAZIONE 

650 INPUT"!M<: (SINISTRA) 3 ";XI 

660 INPUT“X (DESTRA) 3 ",X2 

670 IF X1>=X2 THEN €40 

680 INPUT"MV (BASSA) 3 ",VI 

690 INPUT"V (ALTA) 3 ";V2 

700 IF V1>=V2 THEN 680 

710 FRINT"OSI SONO ASINTOTI VERTICALI ?"j 

720 GOSUB 850PRINTG* IF Qfi="N" THEN NA=-1:GOTO 770 

730 FRINT"«GUANTI ?"i ■GOSUB 850 : PRINTQ 

740 Nfl=Q:FOR K=1 TO NA 

750 PRINT"BASCISSA"K"* ASINTOTO =", : INPUT AS(K) 

760 NEXT K 
770 GOTO 180 


780 REM - 

790 REM STAMPA FUNZIONI 

800 REM - 

310 REM - TRACCIAM. CORNICE 


820 FOR K=0 TO 319 STEP 3:V*0:PLOT X,V,1 : V 3 199 PLOT K,V,1UJEXT K 
830 FOR V=0 TO 199 STEP 3:K=0:PLOT K, V, 1 ■’ K=319 : PLOT K.-V, 1 : NEXT V 

840 COPV-OPEN 4,4 

841 PRIHT44 PRINT44■PRINT44,SPC<5)"ESTREMI DI VISUALIZZAZIONE 
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842 PRINT#4 : PRINT#4,SPC(25) "X <SINISTRA>="X1 

843 PRINT#4,SPC(25)"X <DESTRA>="X2 

844 PRINT#4,SPC < 25 >"V <INFERIORE)=“VI 

845 PRINT#4,SPC(25)"V (SUPERIORE)="V2 

846 PRINT#4 CMD4LIST 450-579 : PRINT#4 : END 

850 REM - 

860 REM LOOP DI ATTESA 

870 REM - 

880 GET Q $-IF 0*="" THEN 880 
890 OQ=ASC<0* > ^ Q=VAL <Qt) RETURN 

900 REM - 

910 REM ELABORAZIONE COORDINATE 

920 REM - 

936 IF V<=V1 THEN 970 
940 IF V>=V2 THEN 970 
950 V=200-(V-V1)/DV 
960 PLOT K, V • 1 
970 RETURN 

980 REM - 

990 REM ASSE X 

1000 REM - 

1010 IF SGN(V1> =3GN(V2) THEN RETURN 
1020 V*INT(V2/DV+.5) 

1025 CHAR 305/7+5,24,1,1 
1030 FOR K=0 TO 319 STEP 3 
1046 PLOT K,V,1 
1050 NEXT K RETURN 

1060 REM - 

1070 REM ASSE V 

1080 REM - 

1090 IF SGN < X1> *SGN(X2) THEN RETURN 
1100 K*-INT(X1/DX+.5) 

1105 CHAR K+5,5,25, t >1 
1110 FOR V=0 TO 199 STEP 3 
1120 PLOT K,V,1 
1130 NEXT V RETURN 

1140 REM - 

1150 REM ASINTOTI 

1160 REM - 

1165 IF NA=-1 THEN RETURN 

1170 FOR J=1 TO NA : IF X>AS(J>-.001 AND X<AS<J>+.001 THEN 1190 
1180 NEXT J RETURN 

1190 FOR V=0 TO 199 STEP 6^PL0T K,V,1=NEXT V 
1200 X=X1+(K+1)*DX: RETURN 


L’alta risoluzione si ottiene nella riga 300 con il solo comando: 
HIRES 6,3 

Il primo numero stabilisce il colore della cornice del video mentre il 
secondo quello della parte centrale dello schermo. 

I colori possono essere variati a piacere usando il codice seguente 
leggermente differente da quello del Basic standard. 


0 = nero 
2 = rosso 


1 = bianco 
3 = cyan 
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4 = porpora 

5 = verde 

6 = blu 

7 = giallo 

8 = arancione 

9 = arancio chiaro 

10 = rosa 

11= cyan chiaro 

12 = porpora chiaro 

13 = verde chiaro 

14 = celeste 

15 = giallo chiaro 


II ritorno alla bassa risoluzione si ottiene con il comando NRM nella riga 
370. 

Il tracciamento si ottiene con il comando: 

PLOT a,b,c 

I primi due numeri rappresentano le coordinate del pixel interessato, 
mentre l’ultima fissa le modalità di tracciamento cioè: 

C = 0 il pixel viene spento 
C = 1 il pixel viene acceso 
C = 2 il pixel viene acceso se spento o viceversa 

Le righe contenenti tale comando sono: 820, 838, 160, 1040, 1120, 
1190. 

II trasferimento del grafico dal video sulla stampante si ottiene sempli¬ 
cemente con il comando COPY alla riga 840. 

Nelle righe seguenti vengono scritti gli estremi di visualizzazione e il 
listato delle linee contenenti le funzioni grafica te. 

Relativamente al tracciamento degli asintoti può accadere che pur senza 
verificarsi un errore causato da una divisione per zero il programma si 
interrompa perché l’ordinata diviene talmente grande da generare un 
errore di overflow. 

Questo raro inconveniente può essere il più delle volte evitato imponen¬ 
do il tracciamento dell’asintoto verticale non appena l’ascissa (assoluta) 
della funzione cada entro un opportuno intervallo contenente l’ascissa 
dell’asintoto. 

Più precisamente, chiamando con XA l’ascissa di un asintoto verrà 
forzato il tracciamento non appena risulterà 

XA - .001 < X < XA + .001 

A questo pensa la riga 1170. 

La struttura del resto del programma è molto simile a quella realizzata 
nel primo capitolo. 
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RICERCA DEI MASSIMI E MINIMI DI UNA FUNZIONE 


Accingiamoci ora a realizzare un programma che consenta la ricerca 
automatica dei punti caratteristici di una funzione (massimi, minimi e 
flessi orizzontali) e alla elaborazione di una tabella contenente le ascis¬ 
se, le ordinate, le derivate prime e seconde di un insieme di punti 
discreti contenuti in un certo intervallo. 


ieri rem **#44#*»4444444**44:4444'4*:*i4*44 
ne REM * * 

126 REM * STUDICI DI UNA FUNZIONE « 

136 REM * * 

146 REM *#*#*#44*4**#44f#i#**4:4********* 

156 CLR = PRINT M tCLEflRKBLEU3":PGKE 53281,15 ;POKE 53286,8 

166 PRINT"[CLEAR3DATTI LA FUNZIONE 

176 INPUT "C3 DOWNIV = ",Vf 

186 PRINT"[CLEAR3[3 DOWN3226 DEF FNY(X)«"V* 

185 PRINT"225 YÌ="CHRÌ<34>V*CHRi(34> 

190 PRINT"DOTO 220[HOME3" 

206 FOR K«0 TO 2 -POKE 631+K,13 : NEXT : POKE 198,8 
216 END 

226 DEF FNV (X> *-XT5/5-X 44/4+XT3/3+X12/2-3 
225 V*="-X15/5-XT4/4+X13/3+X t2/2-3" 

236 DEF FNVKX)*<FNV<:X+lE-4)-FNV<X-lE-4>V2E-4 
246 DEF FNV2(X)=(FNV1(X+1E-43-FNV1(X-1E-4)V2E-4 
256 DEF FNAR(X>=INT(X*lE2+.5>/lE2 
266 POKE 198,6 

276 PRINT"[CLEAR3 OK - ORA DIMMI ENTRO QUALE INTERVALLO" 

286 PRINT" VUOI ESAMINARLA " 

296 INPUT "[DOWN3 X (SINISTRAI";XI 
360 INPUT " X <DESTRA>=“;X2 

316 IF XI>=X2 THFN 276 

326 REM - 

336 REM RICERCA MAX MIN * FLEX 

346 REM - 

350 PRINTTAD<9 ') " [ROSSO3 [ 2 DOWN3UN ATTIMO DI PAZIENZAIELEU3 " 

490 DIM Mi < 26 > ; H- 1 : D= SGNtFNVl (XI > > ■S=. 25 ■ XX*X1 

506 IF SGNtFNYl <XX> >OD THEN S=-SX10'D= SGN(FNVKXX)) 

510 IF FNV1(XX)=0 THEN 566 
526 XX=XX+S ; IF XX<X2 THEN 506 
536 GOTO 636■REM USCITA PER FINE 
560 M=FNAP(XX> 

570 IF SGN(FNV2(XX-.1>><> SGN(FNV2(XX+.1>> THEN 590 
575 IF FNV(XX+.1XFNV(XX) THEN 600 
586 IF FNVCXX+.1>>FNV(XX) THEN 616 
585 GOTO 520 

590 W=FNAR(FNV<XX>>: MÌ<H>«="FLESSO = <"+STRi(M>+";"+STRÌ(W)+"> 
595 PRINT :PRINTMi(H>•PRINT:GOTO 626 

600 W=FNAR < FNV < XX )> : Mi ( H > = " MASSIMO - ( "+STR;ì(M)4-" ; "+STRÌ<W>+" 
)" 

665 PR I NT ■' PR I NTMi ( H ) : PR I NT : GOTO 626 

610 W=FNAR(FNV(XX>) : MÌ(H)="MINIMO = < "+STRi<M)4";"4-STRi(W) + ") 
615 PRINT:PRINTMf(H>:PRINT 

620 H=H+1:D= SGNCFNV1(XX+.1)) : S=.25 : GOTO 520 

630 PRINT"[DOWN3 VUOI UN ELENCO CON LE ORDINATE E LE" 

640 PRINT "BERI VATE PRIME E SECONDE 
660 GOSUB 1000 
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670 IF Q*0"S" THEN 900 

680 PF: I NT " [DOWN3VUOI ANCHE LA STAMPA ?":GOSUB 1000 
680 IF OfO"S" THEN 730 
695 GOSUB 2030 

700 REM - 

710 REM ELENCO 

720 REM - 

730 T1="[P0SS03 X FOO F'<X> F"00” 

740 PR I NT " £ CLE AR 3 CDOWN3 " T4 

750 FOR K«1 TO 40 ■'PRINT" £ ARANC3." ; : NEXT 

760 PRINT"£BLEU3£I'0WN3" : IF FL=1 THEN GOSUB 2046 

770 FOR X=X1 TO X2 STEP .25 

780 PRINTX,FNAR(FNV<X)),FNARCFNV1(X)),FNAR<FNV2< X > > 

790 IF FL«1 THEN GOSUB 2050 

600 T*T+1 : IF I NT < TX 16 ) OTX16 THEN 850 

ei0 PRINTTAB(12>"[DOWNK PREMI UN TASTO)" : GOSUB 100O 

820 PR INT " [ CLEAF: 3 £ BONN 3 " T* 

830 FOR K*1 TO 40 PRINT"IARANO3.".'NEXT 
840 PRINT"t BLEU 31DOWN3" 

850 NEXT 

855 IF FL=1 THEN GOSUB 2130 

860 PRINT"CDOWN3ANCORA ?":GOSUB 1000 

870 IF 0*="S" THEN RUN 

900 PRINT"[DOWN30K - FINE LAVORO." END 

1000 REM - 

1010 REM CICLO BI ATTESA 

1020 REM - 

1030 GET C'f lF G*^"" THEN 1030 
1040 RETURN 

2000 REM - 

2010 REM STAMPA 

2020 REM - 

2030 OPEN 4,4 PR INT#4 , CHF'T •' 14). FL* 1 : RETURN 

2040 PF:I NT#4, Tt ; PRINT#4 : FOF’ K=1 TO 40 ; PRINT#4, "# " ; ■ NEXT 

2045 PF:INT#4 PRINT44 : RETURN 

2050 PF:INT#4,X, 

206O PRINT#4,CHRfO6)"20"FNAR<FNV<X)>; 

2070 PRINT# 4,CHRt< 16>"40"FNAR<FNV1 < X>); 

2080 PR 1NT#4, CHRJ (16)" 60 " FHAF: ( F NV2 ( X ) > 

2090 RETURN 

2100 REM - 

2110 REM STAMPA FINALE 

2120 REM - 

2130 PRINT44:FOR K=1 TO 40 :PRINT44,"#"; :NEXT 

2140 PRINT44 : PRINT44 •' PRINT44, SPC< 10) "V ="VT 

2150 PRINT44 : FOR K=J TO H 

2160 PRINT44,Mf OO 

217G NEXT K 

2180 RETURN 


Nella introduzione della funzione appare una nuova riga (la 185). 
Quando il programma riparte dalla riga 220 (dopo i tre RETURN 
automatici) le variabili vengono tutte azzerate e quindi la Y$ battuta 
nella riga 170 viene “dimenticata” rendendo impossibile la stampa su 
carta del testo della funzione dopo il grafico. 

In 185 viene quindi forzato Tinserimento della nuova riga 225 conte¬ 
nente la funzione sotto forma di stringa che non sarà più dimenticata in 
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quanto si trova dopo il punto in cui il programma ricomincia a scorrere 
(riga 220). 

Nelle righe 230 - 240 sono definite due nuove funzioni, le quali non 
sono altro che i rapporti incrementali della funzione proposta. 


y'(x ) = 



y"(x) 


y' (x + 10~ 4 ) - y' (x - IO" 4 ) 
2 x IO" 4 


Essendo l’incremento sufficientemente piccolo, con buona approssi¬ 
mazione possiamo confondere i rapporti incrementali con le derivate 
prima e seconda calcolate nel punto jc. 

Nella riga 250 si ha una terza funzione con il solo scopo di arrotondare i 
risultati ottenuti prendendo in considerazione solo due cifre decimali. 
Per la ricerca dei massimi, minimi e flessi la routine interessata va dalla 
riga 490 alla 620 e analizza un intervallo (XI-X2) scelto precedente- 
mente. 

In questo intervallo facciamo scorrere un punto da XI a X2 con incre¬ 
menti successivi di lunghezza S. Indichiamo con XX tale ascissa variabi¬ 
le e prendiamo in considerazione il segno D della derivata prima calco¬ 
lata in XX. 

Fino a quando tale segno rimane invariato si avrà la funzione sempre 
crescente (se D>0) o decrescente (se D<0). Non appena verrà regi¬ 
strato un cambiamento di segno, avremo in corrispondenza un punto di 
minimo, di massimo o di flesso orizzontale. 

Per determinare con maggiore precisione le ascisse di tale punto, non 
appena D si inverte cambiamo segno all’incremento S rendendolo 
inoltre più piccolo: l’ascissa XX tornerà indietro più “lentamente” e 
comincerà a oscillare avanti e indietro fino a quando non venga deter¬ 
minato esattamente il punto in cui la derivata si annulla. 

Tale situazione viene rilevata nella riga 510 dove viene effettuato un 
salto in 560 per l’arrotondamento e il riconoscimento (nelle linee suc¬ 
cessive) delle caratteristiche del punto (massimo, minimo o flesso). 
Al termine di tale riconoscimento, nella riga 620 si ha un ritorno 
indietro per il proseguimento dell’analisi del rimanente intervallo. 
Tale ricerca prosegue fino a quando XX non risulta minore di X2, e al 
verificarsi di tale circostanza si ha in 530 il salto per uscire dalla routine. 
Il riconoscimento dei punti caratteristici avviene nelle righe 570 - 615 
semplicemente confrontando le ordinate della funzione e i segni della 
derivata seconda. 

Siax 0 l’ascissa di uno di tali punti caratteristici: se il segno della derivata 
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seconda in (x 0 ~ 0.1) è differente dal segno calcolato in (x 0 + 0.1) allora 
la funzione ha in x 0 una inversione di concavità e quindi x 0 rappresenta 
l’ascissa di un flesso orizzontale. 

Se invece 

f(x o) ^f(x o + 0.1) 

allora x 0 è rispettivamente un massimo o un minimo con coordinata 
Xof(*o) 

Il contatore H (riga 620) serve per la memorizzazione dei punti caratte¬ 
ristici via via rilevati nella variabile M$(H) dimensionata nella riga 490. 
Il resto del listato non presenta particolari difficoltà di comprensione in 
quanto elabora semplicemente una tabella con ascisse, ordinate, deri¬ 
vate prime e seconde neH’intervallo da XI a X2 con incremento .25. 
Alla riga 2000 si ha la stampa della tabella sviluppata per la stampante 
MPS / 801, ma facilmente modificabile per la MPS / 802 o altre stampanti 
eliminando il CHR$(16) nelle righe 2060 - 2080 che serve (unitamente 
alle due cifre tra virgolette) alla tabulazione. 

Nelle righe finali si ha la stampa dei punti caratteristici precedentemen¬ 
te calcolati. 


GRAFICO DI UNA FUNZIONE E CALCOLO DI UN'AREA 


Mentre il programma precedente non richiede l’uso del Simon’s Basic, 
esso è necessario per quello che stiamo per esaminare. Una novità è 
costituita dalla possibilità di calcolare un integrale definito e di vedere la 
corrispondente area sul grafico colorarsi gradualmente dopo il calcolo. 
La funzione da graficare viene introdotta col solito criterio di autostam¬ 
pa: vengono fissati gli estremi di visualizzazione, forniti gli eventuali 
asintoti e stabilito se si vuole calcolare un’area. 

In caso affermativo viene attivato un flag in 325 e si fissano gli estremi di 
integrazione Al e AF. 

La struttura del programma è sostanzialmente identica a quella già vista 
nel capitolo secondo. 

L’unica differenza riguarda il calcolo dell’area, ottenuto senza partico¬ 
lari algoritmi (Gauss-Seidel ecc...), e la colorazione sullo schermo della 
regione interessata. 
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ICt0 CLR PRIHT"."»" 

110 PRINT'TBATTI LA FUNZIONE : " 

120 INPUT" KftWr' = " . Vi 

130 PRINT"muffili80 DEF FNV(X)=“V* 

148 PR I NT " 190 V*= " CHRi ( 34 ) V*CHR* < 34 > 

150 PRINT"G0T0 18048" 

160 FOR K=0 TO 3•POKE 631+K.13 : NEXT■POKE 198 ,4 
170 END 

180 DEF FNV(X)®SIN(X) 

190 VÌ="$IN(X)" 

195 DEF FNflR (V)®INT(V#1E2+.5)/IE2 
■200 POKE 198;0 : PRINT"D3K - ORA ", 

210 PRINT"SCEGLI GLI INTERVALLI DI" 

220 PRINT TAB < 5 >"VISUHLIZZAZIONE 

230 INPUT(SINISTRA)®"!XI 

240 INPUT"X (DESTRA)»";X2 

250 IF X1>=X2 THEN 210 

260 INPUT"»»' (BASSA) = VI 

270 INPUT"V (ALTA)®";V2 

230 IF VI>=V2 THEN 210 

290 PRINT"UDÌ SONO ASINTOTI VERTICALI ? 

300 GOSUB 1O80 : PRINT QJ IF Qi="N" THEN NA®-1:GOTO 322 

310 PRINT"WOUANTI ? "; :GOSUB 1080 :PRINT 0$■ NA=Q 

320 PRINT : FOR K*1 TO NA PRINT"ASCISSA"K"" ASINTOTO INPUT 

AS(K) : NEXT K 

322 PR I NT " STv'UO I CALCOLARE UN' AREA ? : GOSUB 1080 PRI NT Q$ 

324 IF 0*0"S" THEN 360 

325 FL=1 REM FLAG PER CALCOLO AREA 

326 INPUT "ASCISSA INIZIALE" ■ AI 

327 INPUT”«RISCISSA FINALE" i AF 

330 REM - 

340 REM INIZIO ELABORAZIONE 

350 REM - 

360 AIRES 6;3 

370 DX* ( X2-X1 >/320 : DV= ( V2-V1 > /200 

380 GOSUB 570 GOSUB 660 

390 FOR X=0 TO 319 

400 X0=X1+X*DX-GOSUB 960 

410 V0=FNV(X0) 

420 IF VOOV1 THEN 460 
430 IF V0>=V2 THEN 460 
440 V=200-< V0-V1>/DV 
458 PLOT X.V, 1 
460 NEXT X 

470 GOSUB 1130 GOSUB 750 

480 PRINT" VUOI CAMBIARE GLI INTERVALLI DI" 

490 PRI NT "VISUALI ZZAZ IONE ?" 

500 GOSUB 1080 

510 IF Q$="S" THEN PRINT"71 “;:GOTO 210 

520 PRINT"WWVUOI CAMBIARE LA FUNZIONE ?" 

530 GOSUB 1080 

540 IF Qi="S" THEN RUN 

550 F'RINT"WJOK - FINE LAVORO." 

560 END 

570 REM - 

580 REM ASSE X 

590 REM - 

600 IF SGN(VI)=SGN<V2) THEN RETURN 
610 V®INT(V2/DV+.5) 

620 CHAR 305. V+5..24.1,1 
630 FOR X=0 TO 319 
640 PLOT X,V,1 
650 NEXT X RETURN 

660 REM - 

670 REM ASSE V 
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680 REM - 

690 IF SGN < X1> =SGN <X2) THEN RETURN 
700 X=-1NT < X1/BX+.5 > 

710 CHAR X+5,5,25,1 , 1 
720 FOR V=0 TO 199 
730 PLOT X/V/l 
740 NEXT V : RETURN 

750 REM - 

760 REM RITORNO IN BASSA RISOLUZIONE 

770 REM - 

780 GOSUB 1030 

790 IF THEN GOSUB 810 

800 NRM PRI NT ".T RETURN 

810 REM - 

820 REM STAMPA 

830 REM - 

840 FOR X=0 TO 319 : PLOT X,0,1:PLOT X/199/1 NEXT 
850 FOR V=0 TO 199 : PLOT 0,V,l:PLOT 319/V/1 : NEXT 
860 COPV 
870 0PEN4/4 

880 PRINT44 : PRINT#4 : PRINT#4,SPCC5)"GRAFICO DELLA FUNZIONE 
890 PRINT#4:PRINT#4,SPCC30)"V = "V* 

900 PRINT#4 : PRINT#4/SPC<5)“ESTREMI DI VISUALIZZAZIONE : M 

910 PRINT44 : PRINT44.SPC<25>"X <SINISTRA>=”X1 

920 PRINT44,SPC<25>“X <DESTRA)="X2 

930 PRINT#4,SPC<25)"V < INFERIORE.') = "Vl 

940 PRI NT#4.. SPC < 25 ) " V < SUPERI ORE > = "V2 

942 PRINT#4PRINT#4,SPC<5>"ESTREMI DI INTEGRAZIONE ■ " 

944 PRINT#4PRINT#4/SPC<25)"X1 ="AI 
946 PRINT#4/SPC(25)"X2 ="AF 
950 CLOSE 4 -RETURN 

960 REM - 

970 REM ASINTOTI VERTICALI 

980 REM - 

990 IF NA=-1 THEN 1070 

1000 FOR J=1 TO NA 

1010 IF X0=AS(J) THEN 1030 

1020 NEXT J=RETURN 

1030 FOR V=0 TO 199 STEP 2 

1040 PLOT X,V/1 

1050 NEXT V 

1060 X=X+1 :X0=X1+X#DX 

1070 RETURN 

1030 REM - 

1090 REM CICLO DI ATTESA 

1100 REM -- 

1110 GET OS - IF QS= " " THEN 1110 
1120 Q=VAL<Q*> RETURN 

1130 REM - 

1140 REM INTEGRAZIONE 

1150 REM - 

1160 IF FL=0 THEN RETURN 

1170 SU=< AF-AI)/200 

1180 FOR K=AI TO AF-SU STEP SU 

1190 I N= IN+ < FNV ( K ■' +FN V < K+SU '> > *SU/2 

1200 NEXT K 

1210 R*=STR*<FNAR<IN)> -RS-" SPIREA="+R* 

1220 TEXT 5/5/R»/1,1/8 
1230 FOR X0=AI TO AF STEP DX 
1240 X=<X0-X1>/DX 
1250 AV=INT < V2/DV+.5> 

1255 VV=200-<FNV<X0)-V1)XDV 
1260 LINE X/AV/X/W/1 
1270 NEXT X@ 

1300 FL=0 RETURN 
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La routine interessata al calcolo dell’area è quella che va dalla riga 1160 
alla fine. 

In 1160 viene forzato un RETURN nel caso in cui sia nullo il flag relativo 
all’opzione di calcolo dell’area. 

L’intervallo compreso fra gli estremi di integrazione Al e AF viene 
ripartito in 200 parti uguali (SU) e quindi in un ciclo FOR...NEXT viene 
realizzato il calcolo dello scaloide che approssima accettabilmente l’a¬ 
rea della regione curvilinea. 



Fig. 4.1 


Grazie alla suddivisione effettuata, possiamo scomporre la regione 
interessata in una serie di piccoli trapezi aventi altezza SU e basi 
corrispondenti a f(x) e/( jc+SU). 

La loro area è 

f(x) + f(x + SU) 

2 


x SU 
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Con il ciclo FOR...NEXT queste aree vengono gradualmente immagaz¬ 
zinate nell’accumulatore IN nella riga 1190. 

In 1210 e 1220 si ha la stampa su video del risultato mentre nelle righe 
successive si ha la colorazione dell’area calcolata tracciando semplice- 
mente una serie di segmenti verticali ad ascisse variabili da Al ad AF e 
con lunghezze comprese tra l’asse x e l’ordinata della funzione. 

Si noti che se tali segmenti escono dal quadro dello schermo il program¬ 
ma va in errore. 


ESEMPI 


Nella tabella seguente viene riportato un esempio di uscita su stampante 
del programma relativo alla ricerca dei massimi e dei minimi. Nel caso in 


X F(X) F'<X> F " " C X > 


—4 

1 30 

— 2 1 6 

186.26 

- 3 . 7-5 

13 1 . 44 

— 123 . * 4-4 

1 6 * 4.66 

— 3 . 5 

92 . 81 

— 136.5 

134.11 

- 3.25 

62 . 25 

— 1 0 * 4.81 

115.48 

— 3 

40 

— 28 

94 . 99 

- 2 . 25 

23 . 33 

— 55.69 

29 . 35 

- 2.5 

11.81 

— 32.5 

65 . 19 

- 2 . 25 

4 . 32 

— 23.06 

50 . 4 1 

-2 

0 

— 1 2 

38 . 09 

- 1.25 

— 1.93 

— 3 . 9-4 

26 . 4 

- 1.5 

—2 .19 

1.5 

16.9 

- 1.25 

- 1.32 

■4 . 69 

8.21 

— 1 

0 

6 

2 . 05 

- . 25 

1.5 

5.81 

— 3 . 26 

— . 5 

2.81 

■ 4.5 

- 2 . 03 

— . 25 

3 . 69 

2.44 

— 9 . 23 

0 

•4 

0 

— 10 . 01 

. 25 

3 . 69 

— 2 . - 4-4 

— 9 . 26 

. 5 

2.81 

— 4.5 

— 2.04 

. 25 

1 . 5 

— 5 . 81 

— 3 . 26 

1 

0 

— 6 

2 . 05 

1.25 

— 1.32 

— 4.69 

8.21 

1 . 5 

— 2 .19 

— 1.5 

16.9 

1.25 

— 1.93 

3 . 9*4 

26 . 4 

2 

0 

1 2 

38 . 09 

2.25 

- 4.32 

23 . 06 

50 . 4 1 

2 . 5 

11.81 

32 . 5 

65 . 24 

2.25 

23.38 

55.69 

29 . 35 

3 

40 

28 

94 . 99 

3 . 25 

62 . 25 

10 - 4.8 1 

115.48 

3 . 5 

92 . 81 

1 36 . 5 

134 . 11 

3 . 25 

131 . - 4-4 

123 . ■ 4*4 

165.4 

■4 

1 80 

216 

186.26 


***************************************** 

STUDIO DELLA FUNZIONE : 

V -XT4—5*X12+4 


MIMIMO — C — 1 . SB i —2. 23> 
MASSIMO — < 0; 4> 

MI MIMO < 1.53 ;-2.23) 


Fig. 4.2 
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cui si volesse procedere con un incremento più piccolo (per esempio 0.1) 
bisogna inserirlo al posto di 0.25 nelle righe 490 - 620 e 770. 

Se si desidera ottenere tabulati piuttosto lunghi (superiori ad una pagi¬ 
na) è opportuno attivare la paginazione automatica della stampante 
aggiungendo nella riga 2030 un 

PRINT#4, CHR$(147) 

prima del RETURN. 

Il comando di chiusura del canale di stampa è inutile perché se in 870 
viene deciso di procedere ad un nuovo calcolo il comando RUN chiude 
automaticamente tutti i canali compreso quello di stampa. 

Questo programma può essere rapidamente modificato per il Plus 4 e 
per il C16 sostituendo le righe 200 e 260 con le seguenti: 

200 FOR K=0 TO 2:POKE 1319+K,13:NEXT:POKE 239,3 
260 POKE 239,0 

Le quattro figure seguenti rappresentano invece una uscita su stampan- 


rtRErt- 2 

— s 

V 

A A 


X 


GRAFIC0 DELLA FUNZIONE : 


V » SIN(X) 

ESTREMI DI VISUALIZZAZIONE : 

X <SINISTRA)«-3 

X (DESTRA)» 15 

V (INFERIORE)«-2 

V (SUPCRIORE)- 2 

ESTREMI DI INTEGRAZIONE : 

XI ■ 0 

X2 * 3.14 


Fig. 4.3 
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GRAFICO DELLA FUNZIONE : 


V « SJN(X)/X 

ESTREMI DI VISUALIZZAZIONE : 

X <SINISTRA>=-20 

X (DESTRA)* 20 

V (INFERIORE)*-.5 

V (SUPERIORE)* 1.2 

ESTREMI DI INTEGRAZIONE : 

XI * 3.14 
X2 - 6.28 

Fig. 4.4 


te del programma relativo al calcolo degli integrali, con il C64. (La 
modifica per Plus 4 e C16 verrà trattata nel capitolo successivo.) Si noti 
che nel grafico della funzione 

y = sin (x)/x 

la discontinuità per ;t=0 può essere superata segnalando al calcolatore 
un inesistente asintoto verticale che evita l’errore per divisione per zero. 
La funzione successiva presenta l’elaborazione del calcolo seguente 

, , 3,3 - 1 

sen x cos x de =- x-. - 

24 

Gli estremi di integrazione sono stati forniti nella forma approssimata 



XI = .52 
X2 = 1.04 
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Il risultato finale, anch’esso approssimato, è 


AREA = .17 


AREA- .ET 

A, 

V 

i 

V V 

f\ ?K, 

' v v 


GRAFICO DELLA FUNZIONE = 


V - (SIN(X))T2*C0S(X) 


ESTREMI DI VISUALIZZAZIONE : 

X (SINISTRA)—5 

X (DESTRA)» 10 

V (INFERIORE)—.0 

V (SUPERIORE)» .8 

ESTREMI DI INTEGRAZIONE : 

XI » .32 
X2 - 1.04 


Fig. 4.5 
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GRAFICO DELLA FUNZIONE : 


V « (2*X-l)/(2*Xf3) 
ESTREMI DI VISUALIZZAZIONE : 

X (SINISTRA)—2 

X (DESTRA)- 4 

V (INFERIORE)—2 

V (SUPERIORE)- 2 

ESTREMI DI INTEGRAZIONE •' 

XI - .3 

X2 - 4 


Fig. 4.6 




CAPITOLO CINQUE 


Grafici con il Plus 4 e C16 


GRAFICO DI UNA FUNZIONE E CALCOLO DI UNAREA 


Il programma, con alcune modifiche, è identico a quello relativo al C64. 
Sono state modificate le POKE nelle righe 80 e 510 corrispondenti al 
buffer di tastiera e al relativo puntatore, in quanto situate in locazioni di 
memoria differenti. 

Il programma risulta molto più breve (anche perché manca la routine di 
hard-copy) e gira anche sul C16 senza che sia necessario cancellare le 
REM per guadagnare memoria. 

Una differenza molto importante rispetto al listato del C64, dovuta al 
Basic più evoluto, è la capacità di controllare gli errori (in particolare la 
divisione per zero e Yoverflow) senza interrompere lo scorrimento del 
programma. 

Nel caso si verifichi uno di questi errori il programma traccia automati¬ 
camente l’asintoto verticale che quindi non è più necessario dichiarare 
all’inizio. 


10 REM ******************************** 


20 REM * COPYRIGHT IN BASIC * 

21 REM * VIA SEZ2E 22 - LATINA * 

22 REM * TEL. 0773 /487631 * 

23 REM *-* 

24 REM * STUDIO FUNZIONI E CALCOLO * 

25 REM * DI INTEGRALI * 

26 REM * PER C-16 E PLUS-4 * 


30 REM ******************************** 

40 COLOR 0 ,15,2* COLOR 4,7,5:C0L0R 1,8,6 

50 PR INT’iJ3ATTI LA FUNZIONE * - t INPUT " ; Y* 

60 PRINT"*JfcB|500 DEF FNY<X)-"Y* 
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70 PRINT-GOTO 500*" 

00 FOR K-0 TO 2 * POKE 1319+K,13 i NEXT:POKE 239,3 
80 END 

500 DEF FNY(X)*SIN(X> 

510 POKE 239,0:PRINT",JDK“ 

520 PRINT"SCEGLI GLI INTERVALLI DI VISUALIZZAZIONE" 

530 INPUT"**-* (S INI STRA > * " J X 1 

540 INPUT"X (DESTRA) ■ " ;X2 

545 IF X1>*X2 THEN 530 

550 INPUT **Y' (BASSA ) « - ; Y 1 

580 INPUT"Y (ALTA > * " ;Yc 

570 FL«0:PRINT"*^UOI CALCOLARE UN'AREA ?• 

580 GETKEY A* 

590 IF A<< >"S" THEN 1000 

600 FL-1 : INPUT "ASCISSA INIZIALE*"; A I 

610 1NPUT“ASCISSA FINALE** ì AF 

1000 REM - 

1010 REM VISUALIZZAZIONE 

1020 REM - 

1025 PRINT'a" 

1030 DX*(X2-XI>/320 * DY*(Y2-Y1)/200 * GRAPH1C 1,1 

1040 REM - ASSI CARTESIANI - 

1050 IF SGN<Y1>*SGN(Y2) THEN 1070 

1060 Y* INT(Y2/DY♦.5 >:DRAW 1,0,Y TO 319,Y 

1065 CHAR 1,38,l+Y/8,*X" 

1070 IF SGN(X1> *SGN(X2> THEN 1090 

1080 X* -1 NT ( X 1 ,/DX ♦ . 5 ) » DRAW 1,X,0 TO X,199 

1085 CHAR 1,1 *XSB ,1,"Y" 

1090 FOR K«0 TO 319 
1100 X«X1♦K *DX:TRAP 2700 
1110 Y«FNY(X)« GOSUB 4000 
1 120 NEXT 

1130 IF FL*0 THEN 2000 

1 140 REM- 

1 150 REM I NTEGRALE 

1 160 REM- 

1 170 SU»(AF-AI>/200 
1180 FOR K «A I TO AF-SU STEP SU 
1190 IN=IN*(FNY(K)♦FNY(K +SU >)#SU/2 
1200 NEXT K 

1210 IN* INT ( IN* 100* . 5 >/100 
1220 R*»STR<(IN):R*«“AREA»-fR« 

1230 CHAR 1,1,1 ,R* 

1240 FOR X0 «A I TO AF STEP DX 
1250 X«(X0-X1 )./DX 
1260 A Y* I NT ( Y£/DY♦ . 5 > 

1270 YY»200-(FNY(X0)-Y1)/DY 
1280 ORALI 1 ,X ,AY TO X,YY 
1290 NEXT X0 

2000 GETKEY A*:GRAPHIC 0 

2010 PRINT "ANCORA CON LA STESSA FUNZIONE ? * « GETKEY A* 
2020 IF A*«"S" THEN RUN 500 

2030 PRINT"*:ON UN* ALTRA FUNZIONE ?"* GETKEY A* 

2040 IF A*«"S" THEN RUN 

2050 PRXNT-aDK - FINE LAVORO.* 

2060 END 

2500 REM - 

2510 REM TRACCIAM.ASINTOTO VERTICALE 

2520 REM - 

2700 FOR J-0 TO 199 STEP 6 :DRAW 1,K,J*NEXT 
2710 TRAP!RÉSUMÉ NEXT 

4000 REM - 

4010 REM PLOTTAMENTO 

4020 REM - 

4050 IF Y< *Y1 THEN 4 100 
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4060 IF Y>»Y£ THEN 4100 
4070 Y-200-<Y-Yl>/DY 
4080 DRAUI 1,K,Y 
4100 RETURN 


Scorriamo rapidamente il listato. Nella riga 40 vengono fissati i colori (e 
la luminanza) dello schermo, della cornice e dei caratteri. 

In 500 il nuovo comando GET KEY A$ corrisponde alla linea 

580 GET A$:IF A$-""THEN 580 


che genera un ciclo di attesa. 

In 1030 viene attivata l’alta risoluzione con 

GRAPHIC 1,1 


mentre in 1060 e 1080 il comando 
DRAW 1,a,b TO c,d 

traccia una linea fra i punti di coordinate (a,b) e (c,d) che nel nostro caso 
sono gli assi cartesiani. 

Nelle righe 1065 e 1085 vengono stampati i caratteri X e Y in corrispon¬ 
denza degli assi. 

Il comando che permette al calcolatore di tracciare gli asintoti verticali 
in corrispondenza delle divisioni per zero e degli overflow è il TRAP 
(riga 1100) che invia alla subroutine 2700. 

Per il resto il programma è identico a quello già visto per il C64. 
L’uscita su stampante MPS-802 può essere realizzata con una routine 
presa da una analoga di Danilo Toma pubblicata su Computer Club. 
Per il C16 e il Plus 4 occorre modificare le righe 2000-2010 nel modo 
seguente 

2000 GETKEY A*!IF A*<>-e* THEN 2010 

2005 ORAFI 1,0,0 TO 318,01 ORAFI 1,0,189 TO 319,199 

2007 ORAFI 1,0,0 TO 0,199tORAFJ 1,319,0 TO 319,199 

2009 GOSUB 5000 

2010 GRAPHIC 0IPRINT-ÌJ4NCORA CON LA STESSA FUNZIONE ? * i GETKEY A* 

Inoltre si deve aggiungere la subroutine 


5000 REM - 

5010 REM ROUTIFC DI STAMPA PER MPS-802 
5020 REM - 

5030 DATA 189,0,133,253,199,32,133,254,96,160,7,32,108,3,177 
5040 DATA 253,153,124,3,136,16,248,32,116,3,160,7,162,7,94 
5050 DATA 124,3,106,202,16,249,153,132,3,136,16,241,24,165 
5060 DATA 253,105,8,133,253,165,254,105.0,133,254,86,120 
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3070 DAT* 169,252,37,1,133,1,96,169,3,5,1,133,1,68,96 
3060 RESTORE * TD*0 * FOR K-620 TO 691 
5090 READ N«TD*TD+N«POKE K,N 
5100 NEXT 

5110 IF TD 07874 THEN END «REM ERRORE NEI DATA 
5120 C*«" -«REM 20 SPAZI 

3130 OPEN 4,4 « OPEN 5,4,5*0PEN 6,4,6 
5140 SYS 820 

5150 FOR K-0 TO 24 «PRINT46,CHR*<0> 

5160 B4*C*«F0R J-0 TO 39 
5170 SYS 829 

3160 A4* *"* FOR T*900 TO 907 «A4«A4*CHR4<PEEK<T>)*NEXT T 
3190 PRINT#5,A* 

3200 PRINT#4,B4J*B4«B4*CHR»<32> 

5210 PRINT#4,CHR*<254)«NEXT J 
5220 PRINT#S,CHR4<20)* PRINT44 * NEXT K 
5236 CLOSE 4 « CLOSE 3«CL0SE 6 
5240 OPEN 4,4,10*PRINT#4«CLOSE 4 
5250 RETURN 


La stessa routine può essere usata per la stessa stampante e il C64, con 
l’unica modifica di un DATA: il sesto DATA della riga 5030 deve essere 
trasformato da 32 in 224. 


FUNZIONI POLIVOCHE O SOVRAPPOSTE CON IL Ciò E IL PLUS 4 


Anche questo programma è la trasposizione esatta del programma già 
visto per il C64. Le uniche modifiche riguardano i nuovi comandi per 
l’alta risoluzione e il riconoscimento automatico degli asintoti verticali. 


le REM ******************************** 

20 REM * COPYRIGHT IN BASIC * 

21 REM » VIA SEZ2E 22 - LATINA * 

22 REM * TEL. 0773 ."187631 * 

23 REM *-* 

24 REM « STAILA FUNZIONI POLIVOCHE * 

26 REM * PER C-16 E PLUS-4 » 

30 REM ******************************** 

40 COLOR 0,15,2>COLOR 4,7,5>C0L0R 1,6,6 

50 PRINT'JB 1 - VISUALIZZAZIONE FUNZIONI" 

60 PRINTTAB <5)"2 - CAfBIO DELLE FUNZIONI" 

70 PRINTTAB<5>"3 - SCELTA DEGLI INTERVALLI" 

75 PRINTTAB(5>"4 - USCITA PER FINE LAVORO" 

60 GETKEY A* 

85 IF A»-"4" THEN PRINT"*)K - FINE LAVORO. "iEND 
90 IF ASC<A*><49 OR ASC<A*>>51 THEN 80 
100 A«ASC<A8>-46*ON A GOTO 1000,2000,3000 

1000 REM - 

1010 REM VISUALIZZAZIONE 

1020 REM - 

1030 DX* CX2-X1>/320«DY* <Y2-Y1>/200IGRAPHIC 1,1 

1040 REM-ASSI CARTESIANI- 

1050 IF SGN<Y1> »SGN<Y2 > THEN 1070 

1060 Y*INT<Y2/DY +.5)1DRAW 1,0,Y TO 319,Y 

1065 CHAR 1,38,1+Y/8,"X" 

1070 IF SGN<X1)-SGN<X2> THEN 1090 
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1080 X*-INT(X1/DX♦.3)iDRAW 1,X,0 TO X,199 
1085 CHAR 1, 1+X/e, 1,"Y" 

1090 GOSUB 2490 

1100 GETKEY A**GRAPHIC 0 

1110 GOTO 40 

2000 REM - 

2010 REM CAMBIO DELLE FUNZIONI 

2020 REM - 

2030 PRINT-MBH” 

2040 LIST 2500>2650«END 

2490 FOR K *0 TO 319 

2495 X«X1*K*DX*TRAP 2700 

2500 Y»(X-2)/(X*X*X-2>«GOSUB 4000 

2510 * 

2520 Y«(X*X*X+9)/(X*X-2*X-3>xGOSUB 4000 
2530 * 

2540 * 

2550 * 

2560 t 
2570 * 

2580 « 

2590 f 
2600 I 
2610 I 
2620 * 

2630 t 
2640 l 

2650 rem'mhuii^un 40 
2660 NEXT K 
2670 RETURN 

2700 FOR J»0 TO 199 STEP 6 *DRAU l,K,JiNEXT 
2710 TRAP»RÉSUMÉ NEXT 

3000 REM - 

3010 REM SCELTA DEGLI INTERVALLI 

3020 REM - 

3030 PR INT"iJ5CEGL I GLI INTERVALLI DI VISUALIZZAZIONE" 

3040 INPUT"®* <S IN1STRA) "/XI 

3050 INPUT-X (DESTRA >"t X2 

3060 IF X1>«X2 THEN 3030 

3070 INPUT"»T (BASSA)"JY1 

3080 INPUT"Y (ALTA ) " J Y2 

3090 IF Y1>*Y2 THEN 3070 

3300 GOTO 40 

4000 REM - 

4010 REM PLOTTAMENTO 

4020 REM - 

4050 IF Y< *Y1 THEN 4100 
4060 IF Y>*Y2 THEN 4100 
4070 Y«200-(Y-Y1)/DY 
4080 DRAUJ 1,K,Y 
4100 RETURN 


FUNZIONI IN COORDINATE POLARI 


Un punto nel piano può essere individuato con sistemi di riferimento 
diversi da quello cartesiano. Per esempio, basta fissare un punto O e una 
semiretta r uscente da questo punto. La posizione di un altro punto P 
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qualsiasi del piano può essere allora determinata misurando la distanza 
P da 0 e l’angolo che il segmento PO forma con la semiretta r. 

Il segmento r si chiama modulo, mentre l’angolo t si chiama anomalia. 
Immaginando una sovrapposizione fra tale sistema di riferimento e una 
coppia di assi cartesiani con l’origine coincidente con 0 e l’asse r sovrap¬ 
posto a x, si possono elaborare formule di trasformazione da un sistema 
all’altro. 



Fig. 5.2 
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y = r sin(f) 
x - r cos(/) 

Quadrando e sommando si ha: 

y 2 = r 2 (sin(t)) 2 
x 2 = r 2 (cos (t)) 2 

x 2 + y 2 = r 2 [(sin(f)) 2 + (cos(f)) 2 ] 
x 2 + y 2 = r 2 
r = V x 2 + j 2 " 

mentre invece dividendo membro a membro 


ylx = tang(r) 
t = arctang(y/x) 

Talvolta conviene scrivere l’equazione di una curva utilizzando il si¬ 
stema di riferimento polare, perché l’equazione stessa assume una 
forma più semplice. Alcuni esempi: 


r - costante 

t = costante 
r = -8 sin(t) 


r — 4 sin(f) / cos(t) 


r = 30 (l-cos(r)) 
r = l/(sin(z) — 2cos (t)) 


r = 6 tln 
r = 30 sin(3 t) 


equazione di una circonferenza con cen¬ 
tro nel polo 

equazione di una retta con centro nel 
polo 

con le formule di trasformazione si ot¬ 
tiene 

x 2 + y 2 + 8y = 0 

cioè una circonferenza con centro nel 
punto (0; -4) e raggio 4 

equazione di una parabola con centro 
nell’origine e asse coincidente con l’asse 

y 

la curva rappresentata è una cardioide 

equazione di una retta avente in un ri¬ 
ferimento cartesiano la forma 
y = 2x + 1 

spirale di Archimede 

rosa a tre foglie con centro di simmetria 
nell’origine 
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r = 20 / /(cosW — (sin(/)) 2 equazione di un’ellisse avente nel piano 

cartesiano la forma 
x 2 / 400 + y 2 I 100 = 1 


Il programma per visualizzare tali funzioni è il seguente: 


10 REM ********************************* 


80 REM * COPYRIGHT IN BASIC * 
30 REM * VIA SEZZE 28 - LATINA « 
40 REM » TEL. 0773 /487631 * 
50 REM *-* 


60 REM * FUNZIONI IN COORDINATE POLARI » 

70 REM ********************************* 

60 COLOR 0,15/2>COLOR 4,?,5iC0L0R 1,6,6 
90 PRINT'J" 

100 PRINT-LA FUNZIONE DEVE AVERE L'ANGOLO T COME■ 
110 PRINT * VARIAB ILE INO IPENDENTE * 

120 INPUT"MHUR ■")Yt 
130 PRINT-LMM220 R»“Y* 

140 PRINT”RUN 170S" 

150 FOR K-0 TO 8SP0KE 1319+K,13iNEXT!POKE 239,3 
160 END 

170 POKE 239,0!PRINT"JBBFATTORE DI SCALA 7" 

175 PRINT"*10/1 O 1/1 ? <BATT 1 1115 O H25> " 

177 GETKEY A*:IF A*< " 1 “ OR A*>“2" THEN 177 

178 A=10 ! IF A*="2" THEN A=1 
180 GRAPHIC 1 , 1 

185 DPAW 1,0,100 TO 319,100 

167 DRAW 1,160,0 TO 160,189 

130 DO 

200 LOCATE 160,100i TRAP 250 

210 REM-FUNZIONE - 

220 R«30/SQR<COS<T>t2+4*SIN<T)t2> 

230 REM - 

235 IF R<0 THEN 250 

240 LOCATE R*A; 450-G:DRAW 1,RDOT<0>,RDOT<1) 

250 G=G+1:T«G*«/180 

260 R*INT<R*1000 + .5 >/1000 

270 CHAR1,2,2,”R«"+STR*<R>♦" 

280 CHAR1,2,3,"T«"+STR*<G> 

290 LOOP UNTIL G=360 
300 GETKEY A* 

310 GRAPHIC 0 

320 PRINT" J5BANC0RA CON LA STESSA FUNZIONE ? * 

330 GETKEY A* 

340 IFAiF-'S’ THEN RUN 170 

350 PRINT-5CON UNA FUNZIONE DIFFERENTE ?" 

360 GETKEY A* 

370 IFA*=*S“ THEN RUN 

380 PRINT'MOK - FINE LAVORO." 


Il Plus 4 e il C16 hanno la possibilità di controllare e modificare la 
posizione del cursore anche in alta risoluzione. In tal caso prende il 
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nome di pixel cursore e l’unica differenza rispetto al cursore in bassa 
risoluzione è rappresentata dal fatto che non lampeggia. 

È però di notevole importanza avere la possibilità di posizionare il pixel 
cursore nel punto voluto dello schema, prima di iniziare un tracciamen¬ 
to. Il comando relativo è 

LOCATE X,Y 

dove X ed Y sono le coordinate del punto in cui vogliamo indirizzare il 
pixel cursore. 

Vi è la possibilità di usare tale istruzione nella forma 
LOCATE a;b 

dove a rappresenta la lunghezza e b l’angolo (in gradi) dello spostamen¬ 
to desiderato. 

I cicli FOR...NEXT possono essere sostituiti con la coppia di istruzioni 

DO_ 

LOOP UNTIL... 

che significa testualmente 

ESEGUI 

RIPETI FINO A QUANDO... 

Nel listato, l’introduzione della funzione avviene con il solito inserimen¬ 
to automatico nella riga 220. 

Nell’intervallo tra 190 e 290 viene attivato un ciclo DO...LOOP UN- 
TIL... in cui la variabile T (espressa in radianti) e la variabile G (espressa 
in gradi) compiono una rotazione completa. 

Nella riga 200, a ogni ciclo il pixel cursore viene posizionato al centro 
del video. Il TRAP impedisce l’interruzione dell’elaborazione nel caso 
in cui vi sia un errore nella riga successiva. 

In 240 il pixel cursore viene indirizzato nel punto calcolato nella riga 
220 . 

I comandi RDOT(O) e RDOT(1) servono a conoscere la posizione del 
pixel cursore (in coordinate cartesiane relative allo schermo), mentre il 
comando 

DRAW 1, RDOT(O), RDOT(1) 
serve ad accendere il pixel interessato. 
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In 250 si attuano l’incremento della variabile G e la sua trasformazione 
in radianti. 

Si noti come sia stato necessario calcolare il complemento a 450 di ogni 
angolo prima del tracciamento, per compensare il differente metodo di 
determinazione degli angoli fra il sistema americano ed il nostro: in 
quello usato dal calcolatore gli angoli vengono contati a cominciare 
dalle ore 12 e procedono in senso orario. 

Infine in 270-280 si ha la visualizzazione dei valori di R e T calcolati via 
via durante l’elaborazione. 

Si noti che in 235 viene imposto un salto che evita l’accensione del pixel 
nel caso in cui il modulo sia negativo. 


NUMERI COMPLESSI E LORO RAPPRESENTAZIONE NEL PIANO DI GAUSS 


Dato un generico numero complesso 
Z = a + i b 

possiamo rappresentarlo in un sistema di assi cartesiani, considerando il 
coefficiente reale a come ascissa di un punto e il coefficiente immagina¬ 
rio b come ordinata. 

In tal caso il piano prende il nome di piano di Gauss, l’asse delle ascisse 
viene denominato asse reale e l’asse delle ordinate asse immaginario. 
Ogni numero complesso equivale quindi a un punto sul piano di Gauss e 
talvolta conviene far corrispondere a tale numero complesso il vettore 
che parte dall’origine e termina nel punto di coordinate (a, b). 
Considerando il numero complesso come un vettore possiamo applicare 
a esso le considerazioni già svolte nel paragrafo precedente e determi¬ 
nare per ciascuno di essi i corrispondenti valori del modulo e dell’ano¬ 
malia. 

Sui numeri complessi (e sui vettori) è possibile eseguire le quattro 
operazioni aritmetiche e, calcolato il risultato, visualizzarlo sotto forma 
di vettore. 

Con un programma di questo tipo è possibile mostrare l’applicazione 
della regola del parallelogramma sia per la somma sia per la differenza e, 
nel caso di prodotto o quoziente per un numero immaginario, la rota¬ 
zione del vettore risultante. 
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ie rem *************************»****»** 

ee REM * COPYRIGHT IN BASIC * 

30 REM * VIA SEZZE 33 - LATINA * 

40 REM * TEL. 0773 Z48763I * 

50 REM *-* 

60 REM * PIANO DI GAUSS <C-16 PLUS-4) * 

70 REM ********************************* 

80 COLOR 0 ,15,2'COLOR 4,7,5iC0L0R 1,8,6 
30 PRINT'U-IGRAPHIC £,1 

95 REM - 

96 REM TRACCIAMENTO ASSI 

97 REM - 

100 DRAW 1,0,80 TO 319,80 

103 POR K»0 TO 319 STEP 101FOR J-78 TO 8£ 

104 DRAW 1,K,JiNEXT J,K 
110 CHAR 1,33,11,"R" 

130 DRAW 1,160,0 TO 16e,160 

133 FOP K*0 TO 153 STEP 10:FOR J*158 TO 163 

134 DRAW 1,J,K t NEXT J,K 
130 CHAR 1,18,1,* I“ 

140 DRAW 1,0,153 TO 319,159 

150 REM - 

160 REM INGRESSO DATI 

170 REM-- 

175 PR1 NT "SCALA - 1/10 0 1/1 ? HBATTI ai> O CB "» 

177 GETKEY ABtIF AB<>* 1" AND A*<>"3" THEN 177 

178 PRINTASIFM*1lIF AS**3" THEN FM-10 
190 INPUT "PARTE REALE PRIMO VETTORE *' IR 1 
300 INPUT"PARTE IMMAG1NARIA" I I1 

£10 R»RH 1*1 1 iGOSUE 3000 
315 M1*MlA 1*AlGOSUB 3000 

330 INPUT'PARTE REALE SECONDO VETTORE"JR3 
330 INPUT"PARTE IMMAGINARIA" I 13 
340 R-R31I-13lGOSUB 3000 
345 M3*MlA3*AiGOSUB 3000 

350 REM - 

360 REM SCELTA DEL CALCOLO 

370 REM - 

380 PR1NT-BHMB5CEGLI IL TIPO DI OPERAZ IONE* i" 

390 PRINTTAB< 13)"** * /* 

300 GETKEY A* 

310 IF A»-"+* OR A*-*-" OR A*»"*" OR A*-"/“ THEN 330 
330 GOTO 300 

330 IF AB*"*" THEN GOSUB 40001GOTO 500 
340 IF A**"-" THEN GOSUB 50001 GOTO 500 
350 IF A**"*" THEN GOSUB 60001GOTO 500 
360 GOSUB 7000 

370 REM - 

380 REM RISULTANTE 

390 REM - 

500 GOSUB 30001GOSUB 3000 

510 PR1NT-BBBBIL VETTORE RISULTANTE HA i" 

515 R*INT<R*1000*.5)/10001I*INT<I * 1000*.5>/1000 
530 PRINT"» P. REALE* *R;TAB <16 >"P. IMMAGINARIA-■I 
530 PRINT" MODULO*"MITAB <18 >"ANOMALIA*"450“A 
540 GETKEY A* 

550 rn int"—miai irnrn ?■" 

560 GETKEY A* 

570 IF AB*"S* THEN RUN 
580 PRINT-J" 1GRAPHIC 0 
590 PRINT’lfDK - FINE LAVORO1 END 

3000 REM - 

3810 REM CALCOLO MODULO E ANOMALIA 

3030 REM - 

8030 M*SQR < Rt8♦118 >IM* INT<M*100♦.5 >/100 
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6035 IF R-0 AND I >0 THEN 0-90IGOTO 8060 
0037 IF R-0 AND I<0 THEN 0-2701 GOTO 2060 
0040 0-0TN< I/R >*180/+ i0-INT(0 + . 5 ) 

2045 IF R<0 THEN 0-0*180 
0050 IF 0 >360 THEN 0-0-3601 GOTO 2050 
2055 IF 0<0 THEN 0-0*360: GOTO 2055 
2060 0-450-0:RETURN 

3000 REM - 

3010 REM TR0CC1 AhCNTO VETTORE 

3020 REM----- 

3025 DR0W 1,161,79 
3040 LOCATE M»FM> 0 
3050 X1 -RDOT < 0 > : Y 1 -RDOT < I > 

3055 IF X1<0 OR XI >319 THEN 3070 
3057 IF Y1 <0 OR Y1 > 199 THEN 3070 
3060 OR0W 1,161,79 TO XI ,Y1 
3070 RETURN 

4000 REM--- 

4010 REM SOMMA 

4020 REM- 

4030 R-R1*R2*1-I1+12:RETURN 

5000 REM- 

5010 REM SOTTR0Z IONE 

5020 REM---- 

5030 R-R1-R2*I-11- 121RETURN 

6000 REM - 

6010 REM PRODOTTO 

6020 REM--- 

6030 M-M1FM2I0-01+02 

6035 R-R1*R2-11 * 121I-R1FI2+R2*11 

6040 RETURN 

7000 REM--- 

7010 REM OUOZ IENTE 

7020 REM - 

7030 IF M2-0 THEN M-1000IGOTO 7045 
7040 M-M1/M210-01-02 
7045 DD-I2*I2+R2*R2 

7047 R-<R1*R2*I1*I2>/DD:I-<11*R2-R1 » 12>/DD 
7050 RETURN 


Con il comando GRAPHIC 2,1 viene attivata l’alta risoluzione su un 
reticolato di 320x160 pixel, in quanto le ultime cinque righe in basso 
nello schermo rimangono in bassa risoluzione e possono essere sfruttate 
per la visualizzazione dei messaggi necessari durante lo scorrimento del 
programma. 

Dopo il tracciamento degli assi nella routine 150-245 si ha l’introdu¬ 
zione dei coefficienti (reale e immaginario) dei due numeri complessi 
proposti. 

La subroutine 2000 esegue la trasformazione calcolando il modulo e 
l’anomalia dei vettori corrispondenti, mentre la subroutine 3000 traccia 
il vettore stesso. 

Nella routine 250-360 si ha la scelta dell’operazione e lo smistamento 
della subroutine interessata. 

Infine in 370-590 viene visualizzato il vettore risultante. 

Eliminando tutti i REM il programma gira anche sul C16. In tal caso 
occorre cambiare i seguenti numeri di linea 
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2030 va sostituito con 2000 
3025 va sostituito con 3000 
4030 va sostituito con 4000 
5030 va sostituito con 5000 
6030 va sostituito con 6000 
7030 va sostituito con 7000 





Finito di stampare nel settembre 1985 
presso la Lito Velox ■ Trento 
Printed in Italy 




Il libro affronta iproblemi legati alla realizzazione di grafici delle 
funzioni matematiche: spiega linea dopo linea i criteri che guidarlo 
la scrittura di programmi adeguati per visualizzare funzioni 
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Dalle scuole medie all’università, tutti gli studenti troveranno 
in questo libro come ottenere grafici in alta risoluzione (64 000 punti) 
e come stamparli su carta per mezzo di una comune stampante 
Commodore (MPS/801 o MSP/802). 

I programmi girano sia sul Commodore 64 sia su C16 e PLUS4, 
sfruttando il Basic più evoluto di queste macchine. 
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